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L  INTRODUCTION 


Modern  naval  warfare,  especially  undersea  warfare,  depends  heavily 
on  the  vessel  being  acoustically  silent.  A  major  source  of  radiated  noise  is 
the  vibration  of  shipboard  components.  The  reduction  of  these  vibrations  is 
of  utmost  importance  if  a  ship  is  to  accomplish  its  mission.  One  method  of 
vibration  damping  that  shows  promise  in  damping  over  a  broad  spectrum 
of  low  frequency  vibration  is  constrained  viscoelastic  layer  damping.  The 
constrained  viscoelastic  layer  method  uses  the  high  energy  dissipation 
characteristic  of  viscoelastic  materials  during  periodic  motion  of  shear 
deformation  to  absorb  and  dissipate  the  vibrational  energy  of  the  system  in 
question.  Unfortunately  the  design  and  analysis  of  such  constrained 
viscoelastic  layer  systems  is  difficult,  due  in  part  to  the  following: 

•  The  material  properties  of  viscoelastic  damping  materials  vary 
greatly  with  temperature  and  frequency. 

•  Closed  form  solutions  to  the  equations  of  motion  for  constrained 
layer  system  exist  only  for  beams  and  plates  with  simple 
boundary  conditions. 

•  The  exact  complex  valued  eigenvalue  analysis  for  constrained 
viscoelastic  layer  damping  systems  using  the  finite  element 
method  requires  large  amounts  of  computer  storage  and  CPU  time. 


Johnson  and  Kienholz  developed  the  Modal  Strain  Energy  (MSE) 
method  which  uses  the  ratio  of  strain  energy  for  each  mode  shape  to 
approximate  the  modal  damping  of  a  structure  for  a  given  constrained 
viscoelastic  damping  system  [Ref.  1].  This  method  is  very  attractive  due  to 
its  simple  concept  and  very  useful  because  it  can  be  applied  to  any  general 
cases  with  arbitrary  shape  by  using  the  finite  element  method.  However, 
its  effectiveness  compared  with  experiments  were  reported  for  only  a  few 
cases.  Maurer  examined  the  effectiveness  of  the  MSE  method  for  two 
damped  plate  configurations:  1)  a  simple  sandwich  configuration,  and  2) 
a  plate  with  a  milled  pocket  with  damping  material  inserted  and  a  welded 
cover  plate  acting  as  a  constraining  layer  in  a  previous  work  [Ref.2]. 
However,  he  could  not  verify  for  the  milled  pocket  plate  case  since  the  cover 
plate  warped  and  delaminated  itself  from  the  damping  material  during 
welding,  resulting  in  negligible  damping  [Ref.  2]. 

For  certain  naval  applications  the  components  to  be  damped  will  be 
thick  in  construction  and  may  be  exposed  to  an  unfriendly  environment. 
Therefore,  in  this  research  the  pocket  plate  configuration  and  a  second 
milled  plate  using  a  "floating  element"  in  conjunction  with  constrained 
layer  damping  are  investigated.  These  damping  treatments  are  compared 
with  the  simple  sandwich  type  treatments  consisting  of  single  and  double 
constraining  layers.  Simple  plate  geometries  were  used  to  facilitate  the 
experimental  and  computational  effort.  Therefore,  this  paper  addresses 
the  experimental  testing  and  analysis  of  four  thick  aluminum  plates,  each 
with  a  different  constrained  viscoelastic  layer  damping  treatment.      The 


effectiveness  of  the  MSE  method  is  investigated  for  each  different  damping 
treatment  by  evaluating  its  accuracy  in  modal  damping  value  prediction 
compared  with  experimental  results  and  usefulness  as  a  possible  design 
tool. 


II.  THEORY 

A.       VISCOELASTIC  MATERIAL 

Viscoelastic  materials  of  interest  for  general  naval 
applications  are  polymeric  compounds  made  up  of  long  molecular  chains. 
These  molecular  chains  can  be  strongly,  or  weakly,  linked  together, 
depending  on  their  chemical  composition  and  processing.  The  damping 
characteristics  of  viscoelastic  arise  from  the  deformation  and  recovery  of 
the  polymer  network.  Material  properties  of  a  viscoelastic  material  vary 
with  temperature  and  frequency.  As  such,  the  damping  characteristics  of 
a  system  will  vary  as  its  operating  environment  changes.  [Ref.  3] 

Temperature  will  have  the  greatest  effect  on  the  material  properties 
of  damping  materials  [Ref.  3].  This  effect  is  shown  in  Figure  2.1,  where 
four  distinct  regions  are  observed.  The  lowest  temperature  region  is  the 
glassy  region  where  the  material's  storage  modulus  is  at  its  maximum 
value,  and  the  loss  factor  is  at  a  minimum.  In  the  glassy  region  the 
modulus  decreases  slowly  with  temperature  increase,  whereas  the  loss 
factor  increases  rapidly  with  temperature.  The  second  region  is  the 
transition  region  where  the  modulus  decreases  rapidly  with  increasing 
temperature  and  the  loss  factor  reaches  its  peak  value.  The  third  region  is 
the  rubbery  region  where  both  the  modulus  and  loss  factor  are  at  low  values 
and  show  little  variation  with  temperature.    The  fourth,  and  last,  region  is 


the  flow  region  and  characterizes  the  behavior  of  some  materials,  mostly 
ceramics,  at  high  temperatures.  It  should  be  noted  that  the  transition 
region  may  vary  in  width  from  20  °C  up  to  a  width  of  200  °C.  [Ref.  3] 

The  effect  of  frequency  on  viscoelastic  materials  is  not  as  great  as  that 
of  temperature.  The  modulus  of  the  viscoelastic  always  increases  with 
increasing  frequency.  The  loss  factor  will  initially  increase  with  frequency, 
then  peak,  and  subsequently  decrease  as  frequency  increases.  A  plot  of 
storage  modulus  and  loss  factor  versus  frequency  is  shown  in  Figure  2.2.  It 
should  be  noted  that  this  plot  is  over  a  range  of  approximately  ten  decades, 
and  hence  it  becomes  obvious  that  a  temperature  change  of  a  couple  degrees 
will  have  a  much  greater  effect  on  damping  than  a  minor  change  in 
frequency.  [Ref.  3] 

Linear  viscoelastic  materials  behave  in  a  hysteretic  manner  under 
cyclic  excitation.  Therefore,  the  mechanical  behavior  of  a  viscoelastic 
material  during  steady  state  vibration  is  best  described  by  using  a  complex 
stiffness,  k*  [Ref.  3]. 

k*  =  k(l  +  in)  (2.1) 

where, 

rj  =  material  loss  factor 

The  use  of  a  complex  stiffness  then  leads  to  the  use  of  a  complex  Young's 
modulus  and  shear  modulus  [Ref.  3]. 


E*  =  E(l  +  iTi)  (2.2) 

G*  =  G(l  +  rn)  (2.3) 

This  concept  of  the  complex  modulus  is  used  in  subsequent  analysis. 

Viscoelastic  material  properties  are  commonly  displayed  using  a 
"reduced  frequency  nomogram."  The  reduced  frequency  nomogram 
displays  the  variation  of  the  viscoelastic  material's  loss  factor  and  modulus 
with  temperature  and  frequency.     The  "reduced  frequency",  fat,  is  an 

empirically  determined  function  that  accounts  for  the  viscoelastic's 
temperature  and  frequency  dependence,  and  allows  data  for  wide  range  of 
temperature  and  frequencies  to  be  plotted  on  the  same  graph  [Ref  4],  The 
reduced  frequency  nomogram  for  3M  ISD  -  112  is  shown  in  Figure  2.3.  To 
find  the  loss  factor  and  modulus  using  the  nomogram,  enter  with  the 
desired  temperature  and  frequency.  Follow  the  frequency  line  horizontally 
and  the  temperature  line  diagonally  down  the  page  until  the  two  intersect. 
Then  go  vertically  up  or  down  to  intersect  the  shear  modulus  or  loss  factor 
curves.  Finally,  read  the  value  of  the  shear  modulus  or  loss  factor 
horizontally  from  the  scale  on  the  left  [Ref.  4]. 


R        CONSTRAINED  VISCOELASTIC  LAYER  DAMPING 

A  simple  constrained  layer  damping  treatment  consists  of  a  base 
layer  (the  structure  to  be  damped),  a  damping  layer,  and  the  constraining 


layer.   This  configuration  is  shown  in  Figure  2.4  with  the  thicknesses  of  the 
damping  and  constraining  layers  exaggerated  for  clarity. 

The  physical  mechanism  of  damping  can  be  explained  by  referring  to 
Figure  2.4.  When  the  base  layer  is  deformed  in  a  mode  of  vibration,  the 
surface  away  from  the  neutral  axis  elongates,  stretching  the  viscoelastic 
material.  The  top  layer,  being  a  stiff  elastic  material,  tends  not  to  elongate, 
and  thereby  "constrains"  the  viscoelastic  material.  Consequently,  the  cyclic 
motions  of  vibration  induce  a  cyclic  shearing  strain  in  the  viscoelastic. 
This  cyclic  shearing  strain,  together  with  its  associated  hysteresis  loop 
cause  the  vibrational  energy  to  be  dissipated  as  heat.  For  the  constraining 
layer  to  be  effective,  its  stiffness  should  not  exceed  that  of  the  base  layer. 
[Ref.  4  &  5] 


C.        SYSTEM  EQUATIONS  OF  MOTION 

Continuous  systems,  such  as  plates,  possess  distributed 
characteristics  of  mass,  damping,  and  stiffness.  Classical  vibration 
analysis  of  such  systems  involves  the  formation  of  a  mathematical  model 
that  discretizes  the  system  into  a  finite  number  of  components  in  order  to 
approximate  the  total  system.  Such  a  formulation  results  in  the  following 
equation: 

[M]  {x(t))  +  [C]  (x(t))  +  [K]  |x(t)}  =  |F(t)}  (2.4) 


where, 

[M]  =  system  mass  matrix 
[C]  =  system  damping  matrix 
[K]  =  system  stiffness  matrix 
(F(t)}  =  external  excitation  vector 
(x(t)}  =  displacement  vector 

For  an  undamped  system  without  excitation,  the  above  equation  reduces  to 
the  eigenvalue  problem. 

[M]{x(t)}+[K]{x(t)}=0  (2.5) 

This   equation   is   then   transformed   to   modal   space   using  the   linear 
transformation: 


(x(t)}=L<J)J{q(t))  (2.6) 

where, 

[<J>]  =  modal  matrix 

(q(t)}  =  modal  response  vector 

Using  this  linear  transformation  the  equation  of  motion  can  then  be  solved 
for  the  undamped  modal  frequencies  and  mode  shapes. 

To  solve  for  the  frequency  response  of  a  damped  system,  the  linear 
transformation  is  applied  to  equation  (2.4): 

[M]  [<|>]  (q(t)}  +  [C]  [<D]  (q(t)}  +  [K]  [(J)]  (q(t)}  =  (F(t)}  (2.7) 
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Assuming  that  the  damping  matrix  [C]  is  proportional  to  a  linear 
combination  of  the  stiffness  matrix  [K]  and  mass  matrix  [M],  the  damping 
matrix  can  then  be  diagonalized  using  the  same  linear  transformation 
used  to  diagonalize  [K]  and  [M]  in  equation  (2.5)  above.  The  diagonal  terms 
of  the  damping  matrix  then  become  (riiWi),  where  rji  equals  the  modal  loss 
factor  and  cdj  is  the  natural  frequency  of  the  i"1  mode  [Ref.  1]  .  Using  this 
approximate  diagonal  damping  matrix  results  in  a  system  of  uncoupled 
modal  equations  of  motion: 


where, 


qi(t)  +  r\mm)  +  cofaCt)  =  fi(t)  (2.8) 

(j(t)  =  modal  acceleration  of  r"  mode 
qj(t)  =  modal  velocity  of  i"1  mode 
cjj(t)  =  modal  displacement  of  i"1  mode 
"Hi  =  modal  loss  factor  of  i"1  mode 
coi  =  i"1  natural  frequency 
fi(t)  =  modal  force  in  itn  mode 

Assuming  that  a  sinusoidal  excitation  produces  a  sinusoidal  response, 

(fU)}  =  (f)eJ<*       {q(t))  =  {Q}ei<*  (2.9) 

the  response  for  the  i"1  mode  is  then  solved  to  be: 


Q* =  2    «» *  • (2-10) 

c^  -  or  +  jco,rii(jL) 
Subsequently,  the  response  of  the  physical  system  can  be  found  using: 

(x(t)>  =  [(|)]{Q}eJMt  (2.11) 


D.        MODAL  STRAIN  ENERGY  METHOD 

The  equations  of  motion  used  to  define  the  response  of  a  system  with 
viscoelastic  materials  need  a  complex  eigenvalue  analysis.  However,  the 
actual  solution  of  these  equations  may  be  quite  difficult.  This  is  especially 
true  when  the  system  to  be  analyzed  is  comprised  of  materials  whose 
properties  vary  with  both  temperature  and  frequency.  Finite  element 
techniques  are  generally  used  to  compute  the  response  of  complicated 
systems.  However,  for  the  case  of  varying  material  properties  many  time 
consuming  and  costly  runs  must  be  made  with  the  material  properties 
changing  at  each  frequency  increment  [Ref.  1].  In  addition  to  the  costly 
analysis  of  a  single  design  configuration,  changes  in  design  options,  design 
requirements,  of  the  search  for  an  optimum  design  can  make  the  expense 
of  finite  element  analysis  too  great.  The  development  of  the  Modal  Strain 
Energy  (MSE)  method  by  Johnson  and  Kienholz,  however,  makes  the  finite 
element  analysis  of  complex  viscoelastically  damped  structures  a  viable 
option  [Ref.  1]. 
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The  modal  strain  energies  can  be  obtained  for  finite  element  analysis,  and 
are  a  standard  output  option  of  the  NASTRAN  finite  element  code  [Ref.  7]. 

The  modal  frequency  response  of  the  structure  is  then  calculated 
using  the  modal  loss  factors  found  in  equation  (2.12).  When  computing  the 
modal  frequency  response  of  a  damped  structure,  the  modal  properties  in 
the  system  matrices  are  assumed  to  be  constant.  However,  viscoelastic 
materials  have  storage  moduli  which  are  frequency  dependent.  To  account 
for  this  frequency  dependence,  Johnson  and  Kienholz  devised  the  following 
correction  factor  to  be  applied  to  the  modal  loss  factors  calculated  in 
equation  (2.12)  [Ref.  8]. 


where, 


n(r)  =  n(r>      /GST 
V     Cj2,ref 


"He    =  corrected  modal  loss  factor  at  the  r^n  mode 
G2(fr)  =  viscoelastic  shear  modulus  at  the  rtn  modal 

frequency 
Gr2,ref  =  reference  viscoelastic  shear  modulus  used  in  the 

frequency  response  calculation 
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Figure  2.1.   Variation  of  Viscoelastic  Material 

Properties  with  Temperature  [Ref.  3]. 
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Figure  2.2.   Variation  of  Viscoelastic  Material 
Properties  with  Frequency  [Ref.  3]. 
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Figure  2.3.   Temperature  Frequency  Nomogram  for  3M  ISD  -  112. 
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Figure  2.4.    Single  Constrained  Layer  Configuration. 
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IIL  DESIGN  OF  DAMPED  PLATES 

A .        GENERAL  SPECIMEN  CONFIGURATIONS 

For  experimental  testing  and  analysis  purposes,  four  different,  yet 
related,  constrained  layer  damping  treatments  were  selected  in  addition  to 
an  undamped  "reference"  plate.  Two  of  the  damping  configurations  were 
simple  sandwich  treatments  consisting  of  one  and  two  viscoelastic  layers 
respectively.  Another  treatment  was  the  "pocket  plate"  which  was 
previously  investigated  by  Maurer  [Ref.  2].  The  pocket  plate  was  made  from 
a  solid  plate  which  was  then  milled  to  accept  damping  material  and  a  cover 
plate.  The  final  damping  treatment  was  a  "floating  element"  configuration 
consisting  of  a  solid  plate  milled  to  accept  a  double  layer  of  damping 
material  and  a  welded  cover  plate.  Section  views  of  these  damping 
configurations  are  shown  in  Figure  3.1. 

The  purpose  of  the  pocket  plate  is  to  protect  the  viscoelastic  material 
from  materials  such  as  oil  and  salt  water,  which  may  harm  the 
viscoelastic.  Since  previous  attempts  at  using  a  welded  cover  plate  were 
unsuccessful  due  to  the  heat  of  welding  causing  a  delamination  between 
the  damping  material  and  cover  plate  [Ref.  21,  it  was  decided  to  move  the 
viscoelastic  material  from  the  welding  point  by  recessing  it  into  a  shallow 
pocket  of  its  own  as  shown  in  Figure  3.1. 
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The  "floating  element"  concept  evolved  because  the  welded  cover  plate 
of  the  pocket  plate  configuration  does  not  produce  the  damping  reaction  that 
a  true  constraining  layer  would  provide.  If  the  cover  plate  is  welded  to  the 
surrounding  structure  it  cannot  deform  in  bending  as  much  as  an 
unwelded  constraining  layer,  thereby  causing  a  reduction  in  the  damping 
capability  of  the  system.  By  using  a  piece  of  metal  in  the  milled  pocket  with 
dimensions  slightly  smaller  than  the  surrounding  pocket,  along  with  two 
layers  of  viscoelastic  and  a  welded  cover  plate,  a  true  constraining  layer 
effect  should  be  obtained.   This  configuration  is  shown  in  Figure  3.1. 

In  order  to  approximate  a  possible  system  to  be  damped,  a  plate 
with  large  dimensions  was  selected.  The  dimensions  of  the  plates  used  for 
the  damping  treatments  are  114.3  cm  (45  in)  in  length  and  38.1  cm  (15  in) 
in  width.  In  addition,  possible  naval  applications  for  this  type  of  damping 
would  probably  consist  of  thick  plate  members.  For  this  reason  it  was  also 
decided  that  thick  plates  would  be  used  for  the  damping  treatments.  All  the 
base  plates  and  constraining  layers  were  made  of  a  standard  6061-T6 
aluminum  alloy.  The  design  and  selection  of  viscoelastic  layer,  base  layer, 
and  constraining  layer  thicknesses  is  discussed  in  the  following  sections. 


R        DESIGN  OF  THE  SINGLE  DAMPING  LAYER  CONFIGURATION 

In   an   attempt   to   approximate    system   loss   factors   and   hence 
determine  viscoelastic  and  constraining  layer  thicknesses  for  maximum 
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damping,  a  method  developed  by  Nashif  [Ref.  9]  based  on  an  analysis  of 
simple  sandwich  damping  systems  by  Ross,  Kerwin,  and  Ungar  [Ref.  10] 
was  used.  The  Ross  -  Kerwin  -  Ungar  (RKU)  equations  are  base  of  the 
analysis  of  the  simple  sandwich  system  shown  in  Figure  3.2. 

To  find  the  loss  factors  of  the  damped  system,  the  flexural  rigidity  of 
the  system  must  first  be  determined.    For  the  above  system,  the  flexural 
rigidity,  EI,  is  written  as  [Ref.  9]: 


EI  =  EiHi     E2H|     EgHi     EHD2 
12  12  12 

+  E2H2(H2i  -  D)2+  E3H3(H3i  -  D)2 

Egll  +  E2M2^h21  -  D)  +  E3H3(H3i  -  D) 


12 


H31-D 
1  +  g 


(3.1) 


where, 


D  = 


E2H2(h21  -  i|l)  +  g(E2H2H21  +  E3H3H3i) 
EiHi  +  ^M^-  +  g(EiHi  +  E2H2  +  E3H3) 


(3.2) 


H3i  =  ^^.  +  H2 


(3.3) 


H21  =  Hi^H2 
11  2 


g  = 


G2 


E3H3H2K2 


E  =  Young's  modulus 
G  =  Shear  modulus 


(3.4) 
(3.5) 
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I  =  second  moment  of  area 
H  =  thickness  of  member 
K  =  wave  number 

Subscripts  refer  to  the  layers  labeled  in  Figure  3.1.  No  subscript  refers  to 
the  composite  system 

For  a  simply  supported  plate  the  wave  numbers  and  modal 
frequencies  are  found  using  [Ref.9]: 


172       /       EHdgc 

^-^iw^m;  (3-6) 


K^=(if)2+(nbLF  (3J) 


where, 


a  =  semi-wave  length  of  the  plate 

b  =  semi-wave  width  of  the  plate 

v  =  Poisson's  ratio  of  the  composite  plate 

p  =  density  of  the  composite  plate 

gc  =  gravitational  constant 

To  introduce  damping  into  the  equations  it  is  necessary  to  use  the 

complex  modulus   concept  expressed  in   Section   II.      Substituting  the 

appropriate  expressions  for  the  complex  shear  and  Young's  modulus  into 
equations  (3.1),  and  assuming  that  damping  in  the  base  layer,  (r\^),  is 

small,  and  that  the  extensional  stiffness  of  the  damping  layer  is  small 
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(since  E2  «  E^  and  E2  «  E3),  the  following  expressions  can  be  arrived  at 
[Ref.9]: 


EH3  =  EiH?  +  E3H§  +  -J2-  (a  -  p  -  aW  (3.8) 

c2  +  dz 


EH3n  =  E3H3TI3  +  — &-r  (a  -  P  -  dW  (3.9) 

c2  +  d"1 


where, 

a  =  gEiHiEsHsHiijcf  1  -  r|2r,3)  +  6[t\2  +  Hs)  +  J[c(n2  +  nsj])        (3.10) 

(3  =  EiHiE2H2H3i  [c  +  drj2  +  j(cti2  -  d)]  (3.11) 


a  =  2gE2H2E3H3H21H31|    f. f  2™3  ~^J  +  *  "**>    )    (3.12) 

(+  j[c(2ri2  +  Ti3  -  T12H3)  -  d(l  -  2r,2ri3  -  ti|) 


c  =  EiHi(l  +  g)  +  gE3H3(l  -  n2r,3)  (3.13) 

d  =  gE1H1Ti2  +  gE3H3(Ti2  +  r|3)  (3.14) 

j  =  VT  (3.15) 

These  equations  were  then  applied  to  estimate  the  loss  factors  of  the  simple 
three— layer  sandwich  plate.  The  equations  can  be  simplified  by  assuming 
that  damping  in  the  constraining  layer  (T13)  is  negligible  [Ref.9]. 

Since  the  boundary  conditions  for  the  plate  used  in  this  research 
(free-free-free-free)  do  not  correspond  to  the  simply  supported  conditions  on 
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which  equation  (3.7)  is  based,  modal  frequencies  for  the  free-free  case  were 
estimated  using  results  from  finite  element  analysis  of  a  free-free  plate  and 
equation  (3.6).  The  natural  frequencies  of  an  undamped  plate  with  the 
dimensions  previously  given,  and  a  thickness  of  1.91  cm  (0.75  in)  were 
found  using  a  normal  mode  extraction  in  NASTRAN.  By  substituting  these 
modal  frequencies  into  equation  (3.6)  and  estimate  of  the  wave  parameter 
for  each  mode,  K^m,  was  obtained.  Then,  by  using  an  iterative  procedure 
outlined  in  Reference  [9],  modal  loss  factors  were  estimated  for  different 
layer  thicknesses  over  a  temperature  range  of  0.0  °C  to  37.8  °C  (30  °F  to 
100  °F). 

The  previous  equations  are  easily  programmed  to  compute  loss 
factors  for  a  wide  variety  of  conditions.  The  variation  of  viscoelastic 
material  properties  with  temperature  and  frequency  was  accounted  for 
using  a  curve-fit  to  the  reduced  frequency  nomogram  developed  by  Drake 
[Ref.  11].  The  material  data  for  the  following  curve-fit  equations  is  from  the 
University  of  Dayton  Research  Institute  [Ref.  12]. 


logio(M)  =  logio(ML)  + 


'^nffH 


i  ,  (FROMf 

1    FR    I. 


(3.16) 


log10(ETA)  =  logio(ETAFROL)  

+  lc[A(SL  +  SH)  +  (SL  -SH.il-  {T7a*)}     (3'17) 


loglo(FR)  =  loglo(F)-(52^-^0)  (3.18) 
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A  = 


logio(FR)  -  logio(FROL)  (3  19) 


where, 


and, 


M  =  viscoelastic  modulus 

ETA  =  viscoelastic  material  loss  factor 

FR  =  reduced  frequency  (Hz) 

F  =  frequency  (Hz) 

T  =  temperature  (°F) 

T0=4O°C(1O4°F) 

FROM  =  2.0xl04  Hz 

MROM  =  4.75xl06  Pa  (688.94  psi) 

n  =  0.275 

ML  =  6.0xl04  Pa  (8.7  psi) 

ETAFROL  =  1.08 

SL  =  0.45 

SH  =  -0.55 

FROL  =  5000  Hz 

C  =  2.5 


In  addition  to  the  above  constants,  a  Poisson's  ratio  of  0.49  and  a  density  of 
0.909  gram  per  cubic  centimeter  (0.035  lbm/in^)  was  used  for  ISD-112 
[Ref.12].  The  following  material  properties  were  used  for  6061-T6 
aluminum  [Ref.13], 
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E  =  70GPa(10xl06psi) 

v  =  0.33 

p  =  2.7  gm/cm3  (0.0968  lbm/in3) 

Using  the  previous  equations  and  material  properties,  a  computer 
program  was  written  to  compute  estimated  modal  loss  factors  and 
frequencies  for  a  variety  of  base  layer,  viscoelastic  layer,  and  constraining 
layer  thicknesses.  A  listing  of  this  program  appears  in  Appendix  A. 
Modal  loss  factors  were  computed  for  base  layer  thicknesses  of  9.53  mm 
(0.375  in)  to  19.05  mm  (0.75  in)  in  3.18  mm  (0.125  in)  increments.  For  each 
base  layer  thickness,  the  viscoelastic  thickness  was  varied  from  0.38  mm 
(0.015  in)  to  1.52  mm  (0.060  in)  in  0.38  mm  increments,  and  the 
constraining  layer  thickness  was  varied  from  1.59  mm  (0.0625  in)  to  6.35 
mm  (0.25  in)  in  1.59  mm  increments.  In  addition,  loss  factors  were  also 
computed  for  a  viscoelastic  thickness  of  0.127  mm  (0.005  in).  From  the 
results  of  the  analysis,  a  carpet  plot  [Ref.  3],  was  made  for  each  of  the  base 
layer  conditions.  The  carpet  plot  reflects,  for  the  the  first  mode,  the 
maximum  loss  factor  and  its  corresponding  temperature  for  each 
viscoelastic  layer/constraining  layer  thickness  configuration.  The  carpet 
plot  for  a  base  layer  thickness  of  12.7  mm  (0.50  in)  is  shown  in  Figure  3.3. 
Based  on  the  carpet  plots  and  a  desire  for  maximum  damping,  as  well  as  a 
system  which  could  be  moved  easily,  a  base  layer  thickness  of  12.7  mm(0.50 
in),  a  viscoelastic  thickness  of  0.38  mm  (0.015  in),  and  a  constraining  layer 
thickness  of  6.35  mm  (0.25  in)  was  selected.  The  total  system  thickness  was 
approximately  19.05  mm  (0.75  in).  This  total  system  thickness  would  be 
maintained  for  all  subsequent  damping  configurations. 
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To  maintain  continuity  between  damping  systems,  the  milled  "pocket 
plate"  was  given  a  viscoelastic  thickness  of  0.38  mm  and  a  cover  plate 
thickness  of  6.35  mm  for  a  total  system  thickness  of  19.05  mm  (0.75  in).  In 
order  to  keep  the  heat  of  welding  away  from  the  viscoelastic  material,  the 
ISD-112  was  recessed  into  a  shallow  pocket  as  indicated  previously  in 
Figure  3.1,  and  as  shown  in  the  pocket  plate  system  arrangement  in 
Figure  3.4.   Detail  drawings  of  the  pocket  plate  are  shown  in  Appendix  C. 


C.        DESIGN  OF  THE  DOUBLE  DAMPING  LAYER 

The  design  of  the  double  layer  damping  system  was  accomplished 
using  a  modification  of  the  RKU  analysis  used  in  the  previous  section.  The 
RKU  equations  are  used  by  working  from  the  top  layer  of  the  damping 
system  down  towards  the  base  layer.  As  shown  in  Figure  3.5,  the  H3' 
constraining  layer  along  with  the  H2'  viscoelastic  layer  are  combined  with 
the  HI'  layer  to  form  a  three-layer  system.  Using  the  RKU  equations,  the 
stiffness  of  this  system  is  computed  and  considered  to  be  the  equivalent 
stiffness  of  the  top  three  layers  of  the  total  constrained  layer  damping 
system.  The  top  three  layers  were  then  considered  as  a  single  layer  with 
the  equivalent  stiffness  previously  calculated,  and  the  RKU  equations  were 
again  applied  to  compute  estimated  modal  loss  factors  for  the  entire  double 
layer  damping  system  [Ref.3]. 
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To  maintain  continuity  among  all  the  damping  configurations,  a 
base  layer  thickness  of  12.77  mm  (0.50  in)  was  chosen,  and  a  total  system 
thickness  of  19.05  mm  (0.75  in)  was  maintained.  A  design  for  high 
damping  was  then  selected  by  computing  modal  loss  factors  for  the 
constraining  and  damping  layer  thickness  combinations  shown  in  Table 
3.1.  The  estimated  modal  loss  factors  for  the  first  mode  of  vibration  in  each 
configuration  are  plotted  as  shown  in  Figure  3.6.  A  listing  of  the  program 
used  to  compute  the  loss  factors  is  in  Appendix  B. 

From  the  data  presented  in  Figure  3.6,  viscoelastic  thickness  of 
0.38  mm  (0.015  in)  and  constraining  layer  thicknesses  of  3.18  mm  (0.125  in) 
were    selected    for    the    double    layer    configuration.    This    particular 
configuration  estimates  high  damping  over  a  wider  temperature  range 
than  the  other  thickness  combinations. 


TABLE  3.1.  THICKNESSES  USED  IN  CALCULATION  OF  DOUBLE 
LAYER  MODAL  LOSS  FACTORS. 


Configuration 

1 

2 

4 

5 

6 

H2  (mm) 

0.38 

0.38 

076 

0.38 

1.14 

HI1  (mm) 

3.18 

2.38 

2.38 

2.38 

2.38 

H2'  (mm) 

0.38 

0.38 

0.76 

1.14 

0.38 

H3'  (mm) 

3.18 

3.18 

2.38 

2.38 

2.38 

The  milled  "floating  element"  plate  uses  the  same  viscoelastic  and 
constraining  layer  thicknesses  as  the  simple  two— layer  configuration.    In  a 
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design  similar  to  that  of  the  pocket  plate,  the  floating  element  and  both 
layers  of  viscoelastic  are  recessed  into  a  milled  opening  as  shown 
previously  in  Figure  3.1  and  further  described  in  the  floating  element 
system  configuration  of  Figure  3.7. 
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Figure  3.1.    Four  Damping  Treatment  Configurations. 
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Figure  3.2.   Elements  of  a  Simple  Sandwich  Damping  System. 
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Figure  3.3.   Carpet  Plot  of  Maximum  Loss  Factors  for 
a  Base  Layer  with  HI  =  12.7  mm. 


28 


Figure  3.4.   Arrangement  of  the  Pocket  Plate  Configuration. 
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Figure  3.5.   Configuration  of  the  Double  Constrained 
Layer  System. 
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Figure  3.7.   Floating  Element  System  Configuration. 
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IV.  EXPERIMENTAL  RESULTS 

A.        TESTING  ARRANGEMENT 

Experimental  testing  was  performed  on  each  of  the  four  damping 
configurations  and  the  undamped  reference  plate.  In  order  to  approximate 
the  free-free-free-free  boundary  condition,  each  plate  was  suspended  from 
the  roof  of  the  testing  chamber  using  elastic  cords  as  shown  in  Figure  4.1. 
All  of  the  tests  were  performed  in  a  temperature  controlled  environmental 
chamber  which  enabled  temperatures  to  be  maintained  within  ±1  °C.  The 
primary  component  and  user  interface  was  the  Hewlett-Packard  (HP) 
3562A  Dynamic  Signal  Analyzer  (DSA).  The  HP-3562A  was  used  to  provide 
a  swept  sine  signal  to  a  vibration  generator,  and  analyzed  the  returning 
data  signals.  The  HP-3562A  was  used  to  compute  the  frequency  response 
and  coherence  over  a  range  of  50  Hz  to  1050  Hz  using  the  discrete  Fourier 
Transform  in  the  swept  sine  mode.  Ten  averages  were  performed  at  each 
data  point  using  a  frequency  resolution  of  625  mHz  per  step.  The  source 
level  output  to  the  vibration  generator  was  set  at  1.5  volts. 

A  schematic  of  the  testing  apparatus  is  shown  in  Figure  4.2.  Swept 
sine  source  signals  were  fed  from  the  output  jack  of  the  HP-3562A  to  a 
Wilcoxon  F3  vibration  generator  via  the  piezoelectric  output  of  a  Wilcoxon 
PA7C  power  amplifier.  The  vibration  generator  was  mounted  73.48  cm 
(28.93  in)  from  one  end,  and  12.7  cm  (5.0  in)  from  the  front  edge  of  each 
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specimen  as  shown  in  Figure  4.3.  An  integral  force  transducer  was 
mounted  in  the  base  of  the  vibration  generator  to  measure  the  force  input  to 
the  plate.  This  force  signal  was  then  fed  to  input  channel  one  of  the  DSA 
via  a  PCB  462-A  charge  amplifier.  Plate  accelerations  were  recorded  at 
various  points  using  a  PCB  303A-03  accelerometer  as  shown  in  Figure  4.3. 
Acceleration  data  was  fed  to  input  channel  two  of  the  DSA  via  a  PCB  482A05 
power  supply.  Frequency  response  and  coherence  data  was  then  recorded 
on  disk  for  further  analysis. 

Temperatures  within  the  testing  chamber  were  maintained  using  a 
NESLAB  RTE-8  refrigerated  circulating  bath  which  pumped  fluid  through 
a  small  heat  exchanger  in  the  testing  chamber  as  shown  in  Figure  4.2.  In 
order  to  accurately  monitor  the  temperature  of  the  plates,  a  small 
thermocouple  was  inserted  in  the  base  of  each  plate. 

B.        TESTING  PROCEDURE 

1 .        Undamped  reference  plate 

An  undamped,  reference,  frequency  response  measurement 
was  made  at  a  temperature  of  15.6  °C  (60  °F)  to  set  a  standard  response  by 
which  to  measure  the  effectiveness  of  the  damping  treatments.     The 
undamped  frequency  response  was  recorded  over  a  frequency  range  of 
50  -  1050  Hz  using  a  resolution  of  625  mHz  per  point  in  the  DSA. 
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2.         Damped  plate  measurements 

Frequency  response  measurements  of  the  damped  plates  were 
made  at  a  temperature  of  15.6  °C  (60  °F)  at  several  nodes  on  the  plates  in 
order  to  capture  the  damped  response  of  as  many  modes  as  possible.  A 
representation  of  these  nodes  is  shown  in  Figure  4.3.  Responses  were 
recorded  over  a  frequency  range  of  50  -  1050  Hz  with  a  resolution  of 
625  mHz  per  point  in  the  DSA.  Zoom  measurements  were  also  made  to 
capture  better  data  for  certain  modes.  Modal  loss  factors  were  then 
estimated  from  the  frequency  response  and  coherence  measurements 
using  a  curve— fitting  technique  described  in  Reference  [14]. 

C.        SINGLE  DAMPING  LAYER  RESULTS 

The  single  damping  layer  treatment  was  tested  at  4.44  °C  (40  °F), 
15.6  °C  (60  °F),  and  26.7  °C  (80  °F)  so  that  the  effects  of  temperature  on  the 
damping  treatment  could  be  determined.  A  plot  of  the  single  damping 
layer  frequency  response  at  15.6  °C  is  shown  in  Figure  4.4.  The  single  layer 
damping  treatment  resulted  in  high  damping  with  modal  loss  factors 
ranging  from  0.223  at  53.4  Hz  to  0.091  at  876.6  Hz.  Due  to  the  coupling  of 
modes,  loss  factors  for  all  modes  were  not  measured.  The  frequency 
response  of  the  damped  plate  is  characterized  by  a  frequency  shift  to  the  left 
and  a  smoothing  of  the  frequency  response  when  compared  to  the 
undamped  reference  plate.  The  single  layer  treatment  was  especially 
effective  at  reducing  the  frequency  response  of  a  mode  cluster  between  650 
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and  950  Hz.  The  frequency  band  of  this  cluster  was  shifted  approximately 
200  Hz  with  the  amplitudes  of  the  responses  of  the  modes  being 
dramatically  reduced.  The  single  layer  treatment  was  also  effective  at 
reducing  the  amplitude  of  the  response  peaks  over  the  entire  spectrum  of 
measurement.  On  average,  the  highest  peaks  of  the  frequency  response  in 
the  undamped  condition  were  reduced  by  25  decibels,  a  reduction  of  17.8 
times. 

The  effect  of  temperature  on  the  damping  was  quite  pronounced  as 
shown  in  Figure  4.5.  As  the  testing  temperature  was  decreased,  the 
viscoelastic  layer  became  stiffer  and  damping  levels  were  increased.  A 
comparative  listing  of  the  loss  factors  at  different  temperature  is  in  Table 
4.1  and  a  plot  of  the  modal  loss  factors  is  shown  in  Figure  4.6.  Figure  4.5 
shows  the  trend  of  increased  damping  with  temperature  decrease,  and  a 
corresponding  shift  of  modal  frequencies  to  the  right  as  the  viscoelastic 
becomes  stiffer.  These  changes  are  especially  discernible  at  the  lower 
frequencies. 
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TABLE  4.1.  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 

SINGLE  LAYER  AT  DIFFERENT  TEMPERATURE. 


4.44  °C  (40  °F) 

15.6  °C  (60  °F) 

26.7  °C  (80  °F) 

f(Hz) 

II 

f(Hz) 

11 

f(Hz) 

n 

64.3 

0.217 

53.7 

0.223 

49.0 

0.117 

102.8 

0.223 

89.1 

0.145 

83.5 

0.089 

157.4 

0.203 

138.9 

0.172 

129.6 

0.082 

208.4 

0.183 

188.1 

0.184 

177.9 

0.072 

327.5 

0.158 

301.6 

0.120 

243.9 

0.062 

440.1 

0.144 

424.8 

0.109 

409.4 

0.053 

474.9 

0.205 

476.8 

0.111 

450.3 

0.066 

653.1 

0.079 

608.1 

0.067 

552.9 

0.068 

768.0 

0.135 

644.2 

0.069 

586.7 

0.052 

852.0 

0.093 

722.5 

0.130 

629.3 

0.042 

919.0 

0.104 

817.6 

0.081 

681.4 

0.096 

953.5 

0.073 

876.6 

0.091 

790.1 

858.2 
978.1 

0.049 
0.056 
0.069 
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D.        DOUBLE  DAMPING  LAYER  RESULTS 

The  double  layer  damping  configuration  was  also  tested  at  4.44  °C, 
15.6  °C,  and  26.7  °C.  The  frequency  response  of  this  configuration  at  15.6  °C 
as  compared  to  the  undamped  reference  plate  is  shown  in  Figure  4.7. 
Damping  in  the  double  layer  configuration  is  also  high,  with  modal  loss 
factors  ranging  form  0.301  at  53.3  Hz  to  0.107  at  832.4  Hz.  Due  to  modal 
coupling  loss  factors  for  all  modes  were  not  measured.  The  frequency 
response  of  the  two-layer  configuration  is  also  characterized  by  a  dramatic 
reduction  in  response  amplitude  and  a  frequency  shift  to  the  left.  The  peak 
undamped  responses  were  reduced  by  an  average  of  27  decibels,  or  a 
reduction  of  22.4  times  from  the  reference  condition. 

The  effect  of  temperature  on  the  double  layer  damping  treatment  is 
shown  in  Figure  4.8.  As  with  the  single  layer  case,  damping  in  the  double 
layer  configuration  increased  with  a  decrease  in  temperature.  This 
configuration  also  shows  the  shift  of  modal  frequencies  to  the  right  as 
temperature  decreases  and  the  viscoelastic  becomes  stiffer.  Modal  loss 
factors  for  the  double  layer  configuration  are  listed  in  Table  4.2,  and  are 
plotted  for  comparison  in  Figure  4.9. 

To  compare  the  effectiveness  of  the  single  layer  and  double  layer 
configurations  their  frequency  responses  at  15.6  °C  are  plotted  in  Figure 
4.10  with  a  plot  comparing  modal  loss  factors  in  Figure  4.11.    The  two 
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responses  are  quite  similar  although  the  double  layer  configuration  does 
show  an  increase  of  approximately  22  percent  in  modal  loss  factor. 


TABLE  4.2.  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 

DOUBLE  LAYER  AT  DIFFERENT  TEMPERATURE 


4.44  °C  (40  °F) 

15.6  °C  (60  °F) 

26.7  °C  (80  °F) 

f(Hz) 

n 

f(Hz) 

11 

f(Hz) 

n 

51.9 

0.273 

55.3 

0.301 

49.7 

0.202 

122.5 

0.224 

87.1 

0.215 

81.9 

0.188 

212.7 

0.187 

142.1 

0.217 

127.6 

0.156 

290.5 

0.197 

190.5 

0.212 

172.4 

0.117 

382.7 

0.174 

297.8 

0.139 

278.2 

0.097 

494.2 

0.144 

366.8 

0.154 

317.7 

0.077 

558.2 

0.169 

419.3 

0.125 

396.1 

0.070 

602.2 

0.198 

441.2 

0.100 

428.5 

0.077 

642.4 

0.157 

618.7 

0.098 

602.0 

0.067 

752.8 

0.168 

680.9 

0.096 

638.4 

0.072 

815.4 

0.159 

715.5 

0.060 

738.0 

0.071 

896.3 

0.156 

832.4 

0.107 

802.6 

0.077 

968.3 

0.130 

846.1 

0.050 

39 


E.        POCKET  PLATE  RESULTS 

The  milled  pocket  plate  was  constructed  as  previously  shown  in 
Figure  3.3.  The  constraining  layer,  or  cover  plate,  was  welded  in  place 
using  tack  welds  in  an  attempt  to  keep  the  damping  material  away  from  the 
heat  of  welding,  and  the  cover  plate  from  warping,  instead  of  using  a 
continuous  weld  bead  as  was  previously  attempted  [Ref.  2].  The  cover  plate 
was  welded  to  the  base  at  the  corners,  at  the  midpoint  of  the  short  side  and 
at  three  equally  spaced  locations  along  the  long  dimension  as  shown  in 
Figure  4.12.  Following  welding  the  plate  was  tested  to  ensure  that  the 
viscoelastic  had  not  been  damaged  by  the  heat  of  welding. 

The  pocket  plate  was  tested  at  15.6  °C  (60  °F)  and  the  frequency 
response  is  shown  in  Figure  4.13.  The  response  indicates  that  the 
viscoelastic  layer  was  not  damaged  by  welding  and  that  good  damping  was 
attained.  Modal  loss  factors  ranged  from  0.067  at  62.1  Hz  to  0.090  at  923  Hz. 
Although  damping  is  good,  it  is  approximately  half  that  of  the  single  layer 
configuration.  One  reason  for  this  is  that  the  viscoelastic  material  does  not 
completely  cover  the  base  structure.  Another  reason  is  the  presence  of  the 
welded  cover  plate.  Due  to  the  welded  conditions  the  cover  plate  cannot 
induce  shear  deformation  in  the  viscoelastic  layer  as  well  as  a  true 
constraining  layer,  and  therefore  produces  less  damping  than  the  single 
layer  configuration.  The  effects  of  the  welded  cover  plate  are  especially  felt 
in  modes  below  300  Hz  where  the  frequency  response  is  quite  peaked.  The 
response  curve  becomes  more  rounded  and  the  effects  of  the  damping  layer 
are  seen  as  frequency  increases. 
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The  modal  loss  factors  for  the  pocket  plate  are  listed  in  Table  4.3  and 
are  plotted  in  Figure  4.14.  Even  though  the  damping  is  less  than  the  single 
layer,  the  plate  is  still  adequately  damped  as  shown  in  the  frequency 
response  plot  in  Figure  4.13.  In  this  configuration  the  modal  loss  factors 
remained  relatively  constant  throughout  the  testing  spectrum.  The 
increase  in  modal  loss  factor  values  above  800  Hz  is  due  primarily  to  modal 
coupling,  and  the  measured  modal  loss  factors  in  this  range  are  not 
reliable. 


TABLE  4.3  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 
POCKET  PLATE  AT  15.6°C. 

fCHz)        n  f(Hz)        n 


62.1 

0.067 

93.4 

0.042 

142.0 

0.056 

194.9 

0.060 

308.9 

0.051 

437.5 

0.041 

495.6 

0.036 

565.3 

0.043 

625.7 

0.047 

643.2 

0.049 

687.0 

0.044 

841.9 

0.056 

891.5 

0.077 

923.0 

0.090 
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F.        FLOATING  ELEMENT  RESULTS 

The  milled  "floating  element"  configuration  was  constructed  as 
previously  shown  in  Figure  3.7.  The  cover  plate  was  welded  in  a  fashion 
similar  to  the  pocket  plate  as  shown  in  Figure  4.12.  The  center 
constraining  layer,  or  floating  element,  was  made  slightly  smaller  than  the 
surrounding  structure  thus  allowing  the  floating  element  to  act  as  a  "true" 
constraining  layer. 

The  frequency  response  of  the  floating  element  configuration  at 
15.6  °C  (60  °F)  is  shown  in  Figure  4.15.  The  floating  element  is  quite 
effective  as  the  response  shows  a  good  reduction  in  peak  modal  response. 
Measured  modal  loss  factors  range  from  0.089  at  66  Hz  to  0.064  at  935  Hz.  A 
listing  of  measured  modal  loss  factors  is  in  Table  4.4  and  are  plotted  in 
Figure  4.16.  As  with  the  previous  cases,  the  frequency  response  of  the 
floating  element  configuration  is  characterized  by  a  frequency  shift  to  the 
left  and  a  smoothing  of  the  response  as  frequency  increases. 

In  a  comparison  of  the  pocket  plate  and  floating  element 
configurations,  the  frequency  responses  are  plotted  in  Figure  4.17.  A 
comparison  of  modal  loss  factors  for  the  two  configurations  is  shown  in 
Figure  4.18.  The  two  frequency  response  plots  are  similar,  however,  the 
frequency  response  of  the  floating  element  configuration  is  more  rounded 
than  that  of  the  pocket  plate.  The  major  difference  between  the  two 
configurations  is  seen  in  Figure  4.18.    Modal  loss  factors  for  the  floating 
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element  show  an  average  increase  of  25  percent  over  those  of  the  pocket 
plate.  Reasons  for  this  increase  are  the  added  constraining  effect  of  the 
floating  element  and  additional  layer  of  damping  material  present. 


TABLE  4.4  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 
FLOATING  ELEMENT  AT  15.6°C. 

f(Hz)  n  f(Hz)  ii 


66.0 

0.089 

104.9 

0.040 

154.0 

0.094 

211.4 

0.058 

272.8 

0.063 

323.5 

0.075 

435.8 

0.090 

508.1 

0.089 

540.6 

0.115 

618.7 

0.121 

656.7 

0.100 

708.6 

0.088 

857.6 

0.089 

935.4 

0.064 

43 


SUPPORTING  BEAM 


TEST  PLATE 


Figure  4.1.    Testing  Configuration  in  Testing  Chamber. 
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Figure  4.4.   Frequency  Response  of  the  Single  Layer 
Configuration  at  15.6  °C. 
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Figure  4.6.   Comparison  of  Modal  Loss  Factors  for  the  Single 
Layer  Configuration  at  Different  Temperatures. 
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Figure  4.7.  Frequency  Response  of  the  Double  Layer  at  15.6  °C. 
[ :    undamped  response  , :    damped  response] 
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Figure  4.8.   Frequency  Response  Comparison  of  the  Double 
Layer  Configuration  at  Different  Temperatures. 
[ :  4.44  °C, :   15.6  °C  , :   26.7  °C] 
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Figure  4.9.   Comparison  of  Modal  Loss  Factors  for  the  Double 
Layer  at  Different  Temperatures. 
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Figure  4.10.   Frequency  Response  of  the  Double  and 
Single  Layer  at  15.6  °C. 
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Figure  4.11.   Modal  Loss  Factors  for  Single  and  double  layer 
Configurations. 
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Figure  4.12.  Location  of  Tack  Welds  on  the  Cover  Plate 
of  the  Pocket  Plate  Configuration. 
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Figure  4.13.  Frequency  Response  of  the  Pocket  Plate 

Configuration  at  15.6  °C. 
[ :   reference  plate  , :   pocket  plate  ] 
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Figure  4.14.  Modal  Loss  Factors  for  the  Pocket  Plate 
Configuration. 
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Figure  4.15.   Frequency  Response  of  the  Floating  Element 
at  15.6  °C. 
[ :   undamped  plate  , :   floating  element  ] 
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Figure  4.16.   Modal  Loss  Factors  for  the  Floating  Element 
Configuration. 
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Figure  4.17.   Comparison  of  Frequency  Responses  for  the 

Pocket  Plate  and  Floating  Element  Configurations. 
[ :   pocket  plate  , :   floating  element  ] 
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Figure  4.18.   Modal  Loss  Factors  for  the  Pocket  Plate 
and  Floating  Element  Configurations. 
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V.  FINITE  ELEMENT  RESULTS 

A .        UNDAMPED  REFERENCE  PLATE 

The  first  step  in  the  finite  element  analysis  procedure  was  to  model 
and  analyze  the  undamped  reference  plate  for  its  modal  frequencies  and 
frequency  response.  The  finite  element  model  was  generated  using 
PATRAN,  a  computer  aided  interactive  graphics  program  developed  by 
PDA  Engineering.  PATRAN  is  widely  used  for  conceptual,  preliminary, 
and  detailed  design  and  analysis  of  complex  systems.  One  of  PATRAN's 
major  advantages  is  in  the  interactive  construction  of  finite  element  models 
for  use  by  MSC/NASTRAN,  and  its  ability  to  display  MSC/NASTRAN 
results  in  an  easily  understood  graphic  format  [Ref.  15]. 

The  reference  plate  was  modeled  using  84  plate  (QUAD4)  elements  as 
shown  in  Figure  5.1.  The  QUAD4  element  is  an  isoparametric  element 
with  four  nodes,  one  at  each  corner  of  the  element  [Ref.  7].  A  normal  mode 
extraction  was  then  performed  in  order  to  compare  numerical  results  with 
experimentally  determined  modal  frequencies.  Once  satisfactory 
agreement  between  the  modal  frequencies  calculated  in  NASTRAN  and 
those  obtained  experimentally  was  attained,  a  direct  frequency  response 
calculation  was  performed  in  NASTRAN.  This  frequency  response  was 
then  used  as  the  reference  response  for  further  finite  element  models 
incorporating  viscoelastic  damping  treatments.     The  excitation  for  the 
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frequency  response  calculation  was  a  sinusoidal  force  with  an  amplitude  of 
1.0  applied  at  node  66.  The  response  point  was  node  58  as  shown  in  Figure 
5.1.  These  two  nodes  correspond  to  the  points  on  the  plates  used  in  the 
experimental  portion  of  the  research  where  the  vibration  generator  and 
accelerometer  were  attached. 


B.        SINGLE  DAMPING  LAYER 

The  modeling  of  the  single  constrained  layer  damping  system  was 
done  using  techniques  described  by  Johnson  and  Kienholz  [Ref.  1].  As 
shown  in  Figure  5.2,  the  viscoelastic  layer  is  modeled  using  solid  (HEXA) 
elements,  while  the  base  layer  and  constraining  layer  were  modeled  using 
QUAD4  elements.  The  HEXA  element  is  a  solid,  isoparametric  element 
having  eight  nodes,  one  at  each  corner  of  the  element  with  three 
translational  degrees  of  freedom  at  each  node  [Ref.  71.  The  use  of  solid 
elements  for  the  viscoelastic  layer  allows  the  strain  energy  due  to  shearing 
to  be  adequately  represented.  Plate  elements  are  used  in  the  base  layer  and 
constraining  layer  because  of  their  ability  to  account  for  stretching  and 
bending  deformations.  The  plate  element  allows  its  nodes  to  be  offset  from 
the  plate's  center  to  the  surface  of  the  plate,  coincident  with  the  corner 
nodes  of  the  solid  viscoelastic  elements  [Ref.  11.  Thus,  the  single 
constrained  layer  system  was  modeled  using  only  two  layers  of  nodes,  a 
simple  process  in  PATRAN.    For  the  single  layer  configuration,  a  model 
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having  84  elements  per  layer,  and  an  element  meshing  scheme  shown  in 
Figure  5.2  was  used. 

Once  the  damped  plate  had  been  modeled,  normal  mode  extractions 
were  made  using  MSC/NASTRAN.  Five  separate  runs  were  conducted 
using  the  material  properties  of  ISD-112  at  50,  200,  500,  800,  and  1000  Hz 
and  at  a  temperature  of  15.6  °C  (60  °F).  In  addition  to  the  modal 
frequencies,  the  strain  energy  in  the  viscoelastic  elements  and  the  entire 
model  were  output  from  NASTRAN. 

Since  the  shear  modulus  of  a  viscoelastic  material  changes  with 
frequency,  it  was  necessary  to  estimate  the  actual  modal  frequencies  of  the 
damped  plate  using  an  interpolation  procedure  outlined  by  Johnson  and 
Kienholz  [Ref.  16].  The  first  step  of  the  interpolation  process  was  to  plot  the 
shear  modulus  of  ISD-112  versus  frequency  from  5  to  1000  Hz.  Then,  for 
the  first  mode,  using  NASTRAN  results  based  on  ISD-112  material 
properties  a  50  Hz,  the  first  modal  frequency  predicted  by  NASTRAN  and 
the  corresponding  shear  modulus  were  plotted.  The  same  was  then  done 
using  the  first  natural  frequency  predicted  by  normal  mode  extraction 
based  on  ISD-112  material  properties  at  frequencies  of  500  and  1000  Hz.  A 
curve  was  then  passed  through  these  three  points.  The  point  where  the 
NASTRAN  modal  frequencies  for  the  first  mode  intersected  the  ISD-112 
shear  modulus  curve  was  taken  to  be  the  interpolated  modal  frequency  of 
the  single  damping  layer  configuration.    A  plot  of  the  intersection  of  these 
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two  curves  for  the  first  and  second  modes  is  shown  in  Figure  5.3.    This 
interpolation  process  was  then  repeated  for  each  mode  through  1000  Hz. 

Once  the  interpolated  modal  frequencies  were  found,  the  modal 
strain  energy  equations  (2.12)  and  (2.13)  were  used  to  compute  modal  loss 
factors  for  the  single  layer  configuration.  A  set  of  modal  loss  factors  was 
computed  based  on  the  modal  strain  energies  computed  using  viscoelastic 
properties  at  reference  frequencies  of  50  ,200,  500,  and  800  Hz.  A  set  of 
composite  modal  loss  factors  for  the  modal  frequencies  near  these  reference 
frequencies  was  then  selected.  The  resulting  modal  loss  factors  are  shown 
in  Table  5.1  and  are  plotted  versus  frequency  in  Figure  5.4.  As  seen  in 
Figure  5.4,  the  modal  strain  energy  method  is  predicting  high  damping  for 
this  configuration  with  an  average  modal  loss  factor  of  0.195. 


TABLE  5.1.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 
SINGLE  LAYER  USING  THE  MODAL  STRAIN 
ENERGY  METHOD 

f(Hz)  li  f(Hz) 
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0.249 

95 

0.181 

153 

0.234 

202 

0.212 

285 

0.229 

322 

0.214 

447 

0.200 

463 

0.108 

483 

0.258 

517 

0.237 

586 

0.217 

632 

0.164 

663 

0.155 

712 

0.180 

827 

0.134 

869 

0.152 

883 

0.115 
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Using  the  set  of  composite  modal  loss  factors,  the  modal  frequency 
response  of  the  damped  plate  was  computed  using  MSC/NASTRAN.  Modal 
damping  was  introduced  to  the  model  using  the  SDAMP  option  in  the  Case 
Control  Deck  and  the  TABDMP1  damping  table  in  the  Bulk  Data  Deck  as 
described  in  the  MSC/NASTRAN  Handbook  for  Dynamic  Analysis  [Ref.  17]. 
Since  NASTRAN  uses  a  linear  interpolation  between  points  in  the  damping 
table  to  describe  the  modal  damping  in  the  model  [Ref.  7],  a  simple  curve  fit 
was  applied  to  the  set  of  composite  modal  loss  factors  as  shown  in 
Figure  5.4.  Points  from  this  curve  fit  were  then  used  in  the  NASTRAN 
damping  table.  To  compute  the  modal  frequency  response,  a  unit  excitation 
force  was  applied  at  the  same  node  as  the  undamped  plate,  and  the  node 
used  for  the   response  was  also  the  same  as  the  undamped  plate. 

The  results  of  the  modal  frequency  response  calculations  are  shown 
in  Figure  5.5.  The  dashed  line  represents  the  undamped  reference  plate, 
and  the  solid  line  represents  the  modal  frequency  response  of  the  single 
layer  configuration.  Material  properties  at  200  Hz  were  used  for  the 
ISD-112  damping  material.  The  first  thirty  modes  were  used  in  the  modal 
summation  for  the  response.  A  listing  of  the  MSC/NASTRAN  data  deck 
used  to  compute  the  modal  frequency  response  is  in  Appendix  D. 

The  modal  loss  factors  estimated  using  the  modal  strain  energy 
method  are  compared  to  those  measured  experimentally  for  the  single  layer 
configuration  in  Figure  5.6.   The  estimated  loss  factors  are  greater  than  the 
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experimentally  determined  loss  factors  throughout  the  spectrum  of 
interest,  and  especially  in  the  middle  frequencies. 

The  modal  frequency  response  of  the  single  layer  configuration  is 
compared  to  the  experimentally  measured  frequency  response  in 
Figure  5.7.     The  comparison  was  accomplished  by  normalizing  both  the 
experimentally  determined  frequency  response  and  the  frequency  response 

computed  in  NASTRAN.    Both  responses  were  normalized  using  a  value  of 

in/sec 
1.0  — jT —  .    The  effects  of  the  greater  loss  factors  estimated  by  the  modal 

strain  energy  method  are  obvious  as  the  level  of  the  predicted  response  is 
lower  than  the  measured  response.  The  shift  in  frequency  between  the  two 
curves  is  due  to  the  finite  element  model  being  inherently  stiffer  than  the 
actual  system.  The  correlation  between  the  two  curves  is  especially  good 
below  250  Hz  as  this  is  where  the  differences  between  estimated  and 
measured  modal  loss  factors  are  the  least. 


C.        DOUBLE  DAMPING  LAYER 

The  modeling  of  the  double  constrained  layer  damping  system  was 
accomplished  as  shown  in  Figure  5.8.  The  double  layer  configuration 
consists  of  a  base  layer  modeled  with  offset  QUAD4  elements,  two 
viscoelastic  layers  consisting  of  HEXA  elements,  and  the  top  constraining 
layer  modeled  with  offset  QUAD4  elements.   The  middle  constraining  layer 
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was  modeled  using  three  layers  of  HEXA  elements  in  order  to  give  this 
layer  the  stiffness  necessary  to  act  as  a  constraining  layer.  The  model  was 
meshed  using  60  elements  in  each  layer  as  shown  in  Figure  5.8. 

Using  this  model,  the  modal  strain  energy  method  was  applied  to 
determine  approximate  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  To  determine  the  loss  factors,  normal  mode 
extractions  were  performed  using  reference  frequencies  of  50,  200,  500,  800, 
and  1000  Hz.  A  composite  set  of  modal  loss  factors  for  the  double  layer 
system  was  then  compiled  based  on  the  estimated  modal  frequency's 
relation  to  the  reference  frequency  used  to  calculate  modal  strain  energies. 
This  composite  set  of  modal  loss  factors  is  listed  in  Table  5.2,  and  is  plotted 
in  Figure  5.19  as  a  comparison  to  the  experimentally  determined  modal 
loss  factors  for  the  double  layer  configuration.  The  estimated  modal  loss 
factors  for  the  double  layer  show  high  damping,  but  they  compare  favorably 
with  those  measured  experimentally. 

The  modal  frequency  response  of  the  double  layer  configuration  was 
computed  in  a  manner  similar  to  the  single  layer  in  that  smoothed  loss 
factor  data  was  used  in  the  MSC/NASTRAN  damping  table.  Likewise,  a 
unit  excitation  force  was  applied,  and  the  first  thirty  modes  were  used  in 
the  modal  summation.  The  results  of  the  modal  frequency  response 
calculation  are  shown  in  Figure  5.10.  The  dashed  line  represents  the 
undamped  response,  and  the  solid  line  represents  the  frequency  response  of 
the  double  layer  configuration. 
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TABLE  5.2.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 
DOUBLE  LAYER  USING  THE  MODAL  STRAIN 
ENERGY  METHOD. 

f(Hz)  n  f(Hz)  n 


59 

0.278 

89.5 

0.219 

140 

0.239 

189 

0.213 

258 

0.214 

300 

0.205 

402 

0.179 

427 

0.176 

456 

0.264 

486 

0.242 

544 

0.209 

578 

0.141 

593 

0.131 

659 

0.177 

754 

0.102 

783 

0.076 

802 

0.144 

956 

0.077 

969 

0.111 

9&4 

0.048 

The  frequency  response  calculated  using  NASTRAN  was  compared 
to  the  experimentally  determined  frequency  response  of  the  double  layer 
configuration  as  shown  in  Figure  5.11.  Comparison  of  the  two  frequency 
response  curves  shows  similarity  in  form,  but  a  much  lower  response  level 
for  the  numerically  determined  response.  Once  again  this  could  be  due  to 
the  higher  damping  predicted  by  the  modal  strain  energy  method  and  the 
inherently  higher  stiffness  of  the  finite  element  model. 
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D.        POCKET  PLATE  RESULTS 

The  pocket  plate  configuration  was  modeled  with  the  same  offset 
plate  elements  and  solid  viscoelastic  elements  as  the  single  layer 
configuration.  However,  the  pocket  plate  required  the  modeling  of  the 
milled  structure  around  the  viscoelastic  material  and  the  welds  between 
the  cover  plate  and  milled  plate.  A  representation  of  the  model  is  shown  in 
Figure  5.12.  The  base  structure,  cover  plate,  and  the  structure  immediately 
around  the  cover  plate  was  modeled  using  offset  plate  elements.  The 
viscoelastic  material  and  the  portion  of  the  structure  immediately  adjacent 
to  it  were  modeled  using  the  solid  HEXA  elements.  Since  the  viscoelastic 
material  and  cover  plate  are  physically  separated  from  the  surrounding 
plate,  except  where  the  viscoelastic  is  adhered  to  the  base  structure,  care 
was  necessary  in  creating  the  finite  element  mesh. 

The  model  was  created  in  PATRAN  using  PATRAN's  node  editing 
and  equivalencing  capabilities  [Ref.  18].  This  allowed  the  generation  of  a 
finite  element  mesh  with  two  nodes  at  the  same  geometric  point  in  space. 
Using  this  node  editing  technique,  a  mesh  was  created  which  allowed  the 
viscoelastic  and  cover  plate  to  vibrate  separately  from  the  surrounding 
structure,  yet  at  the  same  time,  keep  the  number  of  elements  and  nodes  in 
the  model  to  a  minimum.  The  welded  points  on  the  cover  plate  were  also 
modeled  using  node  editing  techniques.  At  weld  points  the  finite  element 
node  on  the  cover  plate  was  equivalenced  with  its  corresponding  node  on  the 
base  structure,  resulting  in  a  single  node  and  a  connection  between  an 
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otherwise  separate  base  structure  and  cover  plate.  At  non-welded  points 
on  the  cover  plate  there  were  two  nodes  at  the  same  geometric  point;  one  to 
represent  the  cover  plate,  and  the  other  to  represent  the  base  structure.  The 
model  was  meshed  using  a  5x11  mesh  resulting  in  40  elements  in  the  cover 
plate,  viscoelastic  and  base  layer  as  shown  in  Figure  5.12. 

Using  this  modeling  scheme,  the  modal  strain  energy  method  was 
employed  to  estimate  the  modal  frequencies  and  loss  factors.  Normal  mode 
extractions  were  made  using  viscoelastic  material  properties  at  50,  200,  500, 
800,  and  1000  Hz.  Using  these  reference  frequencies  a  composite  set  of 
modal  loss  factors  was  obtained.  These  loss  factors  are  listed  in  Table  5.3 
and  are  plotted  versus  frequency  in  Figure  5.13.  The  estimated  loss  factors 
give  good  damping  over  the  spectrum  of  interest  with  an  average  modal  loss 
factor  of  0.075.  As  with  the  previous  cases,  damping  values  used  for  the 
modal  frequency  response  calculation  came  from  a  curve— fit  to  the  set  of 
composite  modal  loss  factors. 

The  modal  frequency  response  of  the  pocket  plate  was  computed 
using  the  first  30  modes  and  viscoelastic  material  properties  at  200  Hz.  The 
resulting  estimated  frequency  response  is  shown  in  Figure  5.14.  The 
response  shows  a  definite  frequency  shift  to  the  left  along  with  good 
damping  of  the  frequency  response  when  compared  to  the  undamped 
response. 


71 


The  estimated  modal  loss  factors  and  modal  frequency  response  for 
the  pocket  plate  were  compared  to  those  measured  experimentally.  The 
loss  factor  comparison  is  shown  is  Figure  5.15  and  the  frequency  response 
comparison  is  shown  in  Figure  5.16.  The  estimated  modal  loss  factors  are 
higher  than  those  measured  experimentally,  however,  the  frequency 
responses  compare  quite  favorably  with  each  other.  The  frequency 
response  curve  calculated  through  finite  element  analysis  has  a  lower 
response  level  and  a  frequency  shift  to  the  right  of  the  measured  frequency 
response.  This  is  expected  due  to  the  increase  in  damping  predicted  by  the 
modal  strain  energy  method  and  by  the  fact  that  the  finite  element  model  is 
inherently  stiffer  than  the  physical  system 


TABLE  5.3.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 
POCKET  PLATE  USING  THE  MODAL  STRAIN 
ENERGY  METHOD. 

f(Hz)  n  f(Hz)  ii 


67 

0.046 

109 

0.034 

162 

0.091 

217.5 

0.063 

295 

0.091 

333 

0.047 

453 

0.07 

477 

0.093 

511 

0.113 

588 

0.104 

625 

0.086 

668 

0.089 

713 

0.077 

831 

0.077 

834 

0.065 

868 

0.051 
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Figure  5.1.   Finite  Element  Model  of  the  Undamped 
Reference  Plate. 
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RESPONSE  MEASUREMENT 

I 
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0  64  CM  (0.25  IN)  CONSTRAINING  LAYER 
OFFSET  0  318  CM  (0.125  IN) 
84  QUAD4  ELEMENTS 


0  038  CM  (0.015  IN)  ISD-1 12  VISCOELASTIC 
84  HEXA  ELEMENTS 


1.27  CM  (0  50  IN)  BASE  LAYER 
OFFSET  -0.64  CM  (-0  25  IN) 
84  QUAD4 ELEMENTS 


DIMENSIONS  IN  CENTIMETERS  (INCHES) 


Figure  5.2.   Finite  Element  Representation  of  the 

Single  Constrained  Layer  Configuration. 
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Figure  5.3.    Interpolation  of  the  First  and  Second  Modal  Frequencies 
for  the  Single  Constrained  Layer  Configuration. 
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Figure  5.4.   Estimated  Modal  Loss  Factors  for  the  Single  Layer 
Configuration  with  the  Curve  Fit  Used  for  the 
MSC/NASTRAN  Damping  Table. 
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Figure  5.5.   Calculated  Modal  Frequency  Response  of  the  Single 
Layer  Configuration  Using  NASTRAN. 
[ :    reference  plate  , :    single  layer] 
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Figure  5.6.    Comparison  of  Estimated  and  Measured  Modal  Loss 
Factors  for  the  Single  Layer  Configuration. 
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Figure  5.7.   Comparison  of  Estimated  and  Measured  Frequency 
Response  for  the  Single  Layer  Configuration. 
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Figure  5.8.  Finite  Element  Representation  of  the  Double  Layer 
Configuration. 
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Figure  5.9.   Modal  Loss  Factors  for  the  Double  Layer  Configuration 
as  Determined  from  the  Modal  Strain  Energy  Method 
and  Determined  Experimentally. 
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Figure  5.10.   Calculated  Modal  Frequency  Response  for  the  Double 
Layer  Configuration  Using  NASTRAN. 
[ :   reference  plate  , :   double  layer] 
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Figure  5.11.   Comparison  of  the  Experimentally  Determined 
and  Numerically  Predicted  Frequency  Responses 
for  the  Double  Layer  Configuration. 
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Figure  5.12.  Finite  Element  Representation  of  the 
Pocket  Plate  Configuration. 
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Figure  5.13.   Estimated  Modal  Loss  Factors 
for  the  Pocket  Plate. 
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Figure  5.14.   Estimated  Frequency  Response 

for  the  Pocket  Plate  Configuration. 
[ :   reference  plate  , :   pocket  plate] 
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Figure  5.15  Comparison  of  Experimentally  Measured  and 
Estimated  Modal  Loss  Factors  for  the  Pocket 
Plate  Configuration. 
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Figure  5.16.    Comparison  of  Experimental  and  Predicted 
Frequency  Responses  for  the  Pocket  Plate 
Configuration. 
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VL  CONCLUSIONS 

Constrained  viscoelastic  layer  damping  is  an  extremely  effective 
method  for  reducing  broadband  vibration.  In  each  of  the  experimental 
cases  the  peak  amplitudes  of  frequency  response  were  reduced  by 
approximately  25  decibels,  a  reduction  of  18  times  below  the  undamped 
reference  plate. 

In  a  comparison  of  experimentally  determined  modal  loss  factors  for 
the  four  treatments,  the  double  layer  configuration  yields  the  largest 
damping  and  the  pocket  plate  the  least  as  shown  in  Figure  6.1.  Of 
particular  note  are  the  performance  of  the  pocket  plate  and  floating  element 
configurations.  Although  they  are  not  ideal  configurations  in  terms  of 
"true"  constrained  viscoelastic  layer  damping,  the  damping  levels  achieved 
are  quite  satisfactory.  As  was  previously  reported  in  Section  IV,  the 
floating  element  configuration  yielded  an  average  increase  of  25  percent 
over  the  modal  loss  factors  of  the  pocket  plate.  It  is  also  noted  that  the 
average  modal  loss  factor  for  the  pocket  plate  is  approximately  50  percent  of 
the  average  modal  loss  factor  of  the  single  layer  configuration.  Similarly, 
the  average  modal  loss  factor  for  the  floating  element  configuration  is 
approximately  50  percent  that  of  the  double  layer  treatment. 

The  modal  strain  energy  method  tends  to  overpredict  the  modal  loss 
factors  by  as  much  as  50  percent  for  these  highly  damped,  thick  plates. 
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Figure  6.2  shows  that,  although  the  modal  strain  energy  method  predicted 
modal  loss  factors  greater  than  those  measured,  the  same  relative 
differences  in  modal  loss  factor  between  damping  configurations  are 
maintained.  The  estimated  modal  loss  factors  for  the  pocket  plate  are 
approximately  60  percent  less  than  those  of  the  single  layer  configuration. 
This  indicates  that  there  is  a  consistency  between  the  damping  values 
predicted  by  the  modal  strain  energy  method  and  those  of  the  actual 
physical  system. 

There  are  several  possible  reasons  for  the  differences  between  the 
experimentally  determined  modal  loss  factors  and  those  estimated  by  the 
modal  strain  energy  method.  The  first  is  that  the  material  properties  of 
ISD-112  as  reported  by  the  3M  Corporation  on  the  reduced  frequency 
nomogram  may  not  be  consistent  with  the  material  properties  actually 
present  in  the  material  used.  Since  the  numerical  analysis  was  based  on 
the  reported  material  properties  this  is  a  possible  source  of  uncertainty  in 
the  results. 

The  second  source  of  uncertainty  is  in  the  adhesion  of  the  ISD-112  to 
the  aluminum  plates.  Although  the  plates  were  clean  when  the  ISD-112 
was  applied,  it  was  noted  that  the  adhesive  qualities  of  the  ISD-112  were  not 
uniform  throughout  the  material.  A  lack  of  adhesion  may  cause  a 
reduction  in  the  damping  capability  of  the  system. 
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Another  source  for  the  difference  between  the  experimental  results 
and  numerical  results  lies  within  the  finite  element  model.  The  number  of 
elements  used  in  the  model  greatly  affects  its  "stiffness."  As  the  number  of 
elements  is  increased  the  model  should  become  less  stiff  and  results  are 
expected  to  approach  those  that  are  measured.  Also,  the  type  of  element 
used  to  model  the  base  layer  and  constraining  layer  may  have  an  effect.  In 
this  research  plate  elements  were  used  to  model  both  the  base  and 
constraining  layers.  It  is  possible  that  one  or  more  layers  of  solid  elements 
may  produce  results  that  agree  better  with  experimental  results. 

One  drawback  to  the  modal  strain  energy  method  in  design  is  the 
large  amount  of  CPU  time  required  for  normal  mode  extraction  and  modal 
frequency  response;  especially  in  complex  structures  with  a  large  number 
of  elements.  Therefore,  although  the  modal  strain  energy  method  is  good 
for  analyzing  a  design,  it  may  face  a  big  difficulty  to  be  used  for  design 
optimization  due  to  the  large  amount  of  CPU  time  required  for  normal 
mode  extraction  and  modal  frequency  response  calculations. 
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Figure  6.1.    Comparison  of  Experimentally  Measured   Modal  Loss 
Factors  for  the  Four  Damping  Configurations. 
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Figure  6.2.  Comparison  of  Estimated  Modal  Loss  Factors  for 
the  Single  Layer,  Double  Layer,  and  Pocket  Plate 
Configurations. 
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VII.  RECOMMENDATIONS 

The  thick  plate  used  in  this  research  is  a  generic  model  of  many 
physical  systems  that  may  see  use  in  naval  application.  There  are  several 
areas  which  deserve  more  research  and  clarification,  including  the 
following: 

•  Model  the  floating  element  configuration  in  finite  elements  and 
check  the  effectiveness  of  the  modal  strain  energy  method  in 
predicting  modal  loss  factors  for  this  configuration. 

•  Investigate  the  relationship  between  mode  shapes  and  damping 
values.     It  seems  some  modal  damping  values  are  very  low  due  to 
their  twisting  mode  shapes. 

•  For  pocket  plate  and  floating  element  configurations,  weld  the  whole 
cover  plate  (continuous  weld)  and  investigate  the  damping 
characteristics   to  compare  the  vibration  reduction  capabilities  of 
the  two  different  milled  plate  damping  treatments. 

•  Investigate  methods  for  improving  the  adhesion  of  ISD-112  to  the 
base  structure  and  constraining  layer. 
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APPENDIX  A 

FORTRAN  PROGRAM  USED  TO  COMPUTE  MODAL  LOSS  FACTORS 
FOR  THE  SINGLE  DAMPING  LAYER  DESIGN 

This  program  used  the  Ross-Kerwin-Ungar  equations  of  Section  III 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  single  layer 
configuration.  Modal  frequencies  for  an  undamped  plate  are  read  from  a 
data  file  and  estimated  modal  frequencies  and  loss  factors  for  various  layer 
thicknesses  are  output  to  another  file.  Material  properties  of  ISD-112  are 
computed  using  University  of  Dayton  data  and  curve-fitting  equations  to  the 
reduced  frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program 
are  pounds,  inches,  and  seconds. 
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ITRITE(10,702)  'VISCOELASTIC  LAYER  THICKNESS,  H2=',H2 

IIRHEC  10,X) 

MRU  EC  10,702)  'CONSTRAINING  LAYER  THICKNESS,  H3=',H3 

FORMATC A35,  F6  .  <i  ) 

HRITEUO.x) 

HRITE( 10,700)  'TEMP' ,  'MODE*.  'FCP','ETA2',  'ETAS',  'G2' 

F0RMAT(A8,  3X,A5,4A15) 

PLATE  MODE  LOOP,  READ  MODAL  FREQUENCY  AND  COMPUTE  HAVE  NUMBER 

HTOT  =  HHH2*H3 

DO  300  1=1,8 
READUl,*)  HP(I) 

CALCULATE  MODAL  FREQUENCY  AND  HAVE  NUMBER  OF  BASE  PLATE 

SUB7  =  SQRT((Elx(HTOTxx3)xGC)/(12.x(l  .-NUlxx2)xRH01«lll  )) 
KQR  =  ITP(  I  )/SUB7 
FP=HP(I )/(2.xPI ) 

TEMPERATURE  LOOP 

T=30.0 

DO  <i0  0  K  =  l,  15 

CALCULATE  PROPERTIES  OF  VISCOELASTIC  FOR  GIVEN  TEMPERATURE  AND 
MODE 


C 
C 
C 

c 

501    FR10=LOG10(FP)-(12.x(T-T0))/(525.*:-T0) 
FR=10. xx(FRlO) 


97 


GQR  r ( 1  .  IAXX2.  )  ) )/2 


AM  I'Rl  0-LOO)  0(  r  ROL  )  )/C 
SUB1-Cx(  (SUSII)  xAKSL -f.il  )*<1 
El  A21U  =  L()G10(  E1FRUL  )  •  SUbl 

VISC0ELAS1 ]C  LOSS  I  ACTOR 

ETA2=10.XX(riA210) 
SUD2-2  .  *LOG10(MRUM/ML  ) 
SUD3=1  .1(1  R()H/FR)KX|| 
fllO  =  lOG10(ML)*  5UB2/SUI13 

VJSCOELASTIC  SHEAR  MODULUS 

G2=10.  xx  CI11  0  ) 
VISCOEIASTIC    YOUNG'S    MODULUS 

E2-G2X2.  X(  1  .HIIU2) 

CALCULATIONS  FOR  DAMPED  PLATE  AT  TEMPERATURE  T  USING  ROSS- 
KCRMIN-UNGAR  EQUATIONS 

G=G2/( C3XHJXH2XKQR) 

H2I  =(HHII2)/2. 

I!51=M2+(II1  HI3)/2. 

Cl  =  Clx|llx(  1  .  <G)+G*E3*H3 

D  =  GxETA2X(Flx.Nl  +  E3X113) 

Al  rilRE  =  Gxeix|llxE3xil3x(H3lxx2.  )x(Cl+D*ETA2) 

ALPHIM  =  GxElx||lxE3xH3x(M31X¥2.)xClxETA2 

DRE  =  E1XII1XE2X  112*11 31X(C1  +  D*ETA2) 

BIM  =  Elx||lxE2x|l2xH3lX(ClxETA2-D) 

SUn8  =  2.XGxE2xH2xE3xM3x|l2lxH3lxc 

DELRE=SUB8x((l.-(ETA2xx2.))+(Dx(2.*ETA2))) 

DFLIM=SUB8X(C1X(2.XETA2)-(DX(1.-(ETA2XK2.)))) 

SUB'i  =  (12./(Clx*2.-»DX*2.))x(ALPHRE-nRE-DELRC) 

EHCUDE=(Elx(HlXK3.  )HE3x(M3xx3.)  +  SUB4 


MODAL  FREQUENCY  OF  DAMPED  PLATE 

DEHS  =  RHOlxHHRH02xH2-fH3xRH03 
SUB5=(EHCUDE*GC)/(12.x(l  . -NU1 x*2 . ) XDENS ) 
IICP  =  KQRxSQRl  (SUB5) 
FCP=HCP/(2.XP1 ) 

COMPARISON  OF  FP  AND  FCP 

IF  (ABSd  .-FP/FCP)  .LE.  0.1)  THEN 

GOTO  500 
ELSE 

FP=FCP 

GOIO  501 
ENDIF 

COMPUTE  SYSTEM  LOSS  FACTOR 

00  SUB6  =  (12./(Clxx2.4nxx2.))x(ALPIIIM-BIM-DELIM) 
El AS=(  1  ./EIICUBE)XSUD6 

PRINT  RESULTS 

t IR 1 1  EC  10,701)    T,I,FCF\E1A2,ETAS,G2 

01  I ORMAT(5X,F7  .  3  ,  2X,  I  2,  3X,  '.El  5  .  '.  ) 

NEXT    TEMPERATURE 

T  =  T  +  5.  0 
00         CONTINUE 

NEXT    MODE 

IIRITE(  10,  x) 
iOO  CONI1NUE 
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REIHNIKUNI  r  =  1  1  ) 

c 

C       NEXT    C0II5TP.AII1IIIG    LAYER    THICKNESS 
C 

H.T>  =  II3<  0  .  0625 
200  CONTINUE 

c 

C       NEXT    VISCOELASUC    LAYER    THICKNESS 
C 

112  =  11210.015 
150         COIIlIliUE 
C 

C       NEXT    BASE    TLATE    THICKNESS 
C 

111  =nnn .  125 
100  CONTINUE 

Ct  OSE<  l)NIT  =  l  0) 

CLOSECUNI T  =  1 1 ) 

END 
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APPENDIX  B 

FORTRAN  PROGRAM  USED  TO  COMPUTE  MODAL  LOSS  FACTORS 
FOR  THE  DOUBLE  DAMPING  LAYER  CONFIGURATION  DESIGN 

This  program  uses  the  Ross-Kerwin-Ungar  equations  of  Chapter  3 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  Modal  frequencies  for  an  undamped  plate  are  read 
from  a  data  file  and  estimated  modal  frequencies  and  loss  factors  for 
various  layer  thickness  combinations  in  the  double  layer  configuration  are 
output  to  another  data  file.  Material  properties  of  ISD-112  are  computed 
using  University  of  Dayton  data  and  curve-fitting  equations  to  the  reduced 
frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program  are 
pounds,  inches,  and  seconds. 
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PROGRAM     1IIULYR 

c 

C  *XXXXXXXXXMXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XXX 
C 

C  IHJ5  PROGRAM  COMPOTES  THE  SYSTLM  1053  TAd  OR  AMD  MODAL  FREQUENCY  01' 

C  A  DOUBLE  CONSTRAINED  LAYER  V  I  5C0EL  AS  I  I  C  AL  L  Y  DAMPED  PLAIE.   THE  LOSS 

c  rAciuKn  are  compuied  tor  a  specific  ri  a  i  e/i>ampiiig  layer 

C  CUIIf  Ji'.URAl  IOH  AND  OVER  A  TEMPERAIUR'C  RANGE  01"  50-100  DEGREES 

C  FAKENHE1T. 

C 

C  THIS  PROGRAM  APPL1C3  TO  A  ERE  E  -  F  R  E  E-  FRE E-  T  R  EE  ELATE  AND  HIE  UNDAMPED 

C  MODAL  IREQUENCIES  ARE  DETERMINED  FROM  rillllE  ELEI1EN1  ANALYSIS  AND 

C  ARE  READ  INTO  THIS  PROGRAM  FROM  FILE  'PL1TRQ'. 

C 

C  THE  VISC0EIA3TIC  MATERIAL  IS  311  I  SD-  112.   VISCOELASTIC  MATERIAL 

C  DATA  IS  TROM  UNIVERSITY  OF  DAYTON  RESEARCH  1 1131  I  T  U  T  E 

C 

C  THE  UNITS  U3ZD  III  THIS  PROGRAM  ARE  LD,  INCH,  SEC,  AND  DEGREES  T. 

C 

C  LOSS  FACTORS,  DAMPED  PLATE  MODAL  FREQUENCIES,  AND  I SD- 112  PROPERTIES 

C  ARE  COMPUTED  III  SUBROUTINE  'RKU*  FOR  EACH  ll-TH  C0HSTRA1IIED  LAYER 

C  SYSTEM.   UPON  COMPLETION  OF  EACH  TEMPERATURE  COMPUTATION  THE  SYSTEM 

C  LOSS  FACTOR  AND  CORRESPONDING  DAMPED  PLAIE  FREQUENCY  ARE  MRITTEN  10 

C  TILE  '  1II0LYR  DATA'  . 

C 

C  THE  roi  LOIIIIIG  COEFFICIENTS  ARE  DEFINED: 

C  El  =  YOUNG'S  MODULUS  OF  BASE  PLAIE  (PSI) 

C  E3  =  EQUIVALENT  YOUNG'S  MUDUIUS  OF  N-TH  CONSTRAINED  LAYER  SYS1EM 

C  E3PRM  =  YOUNG'S  MODULUS  OF  THE  CONSTRAINING  LAYER  IN  N-TH  LAYER 

C  EHCUDE  =  EQUIVALENT  STIFFIIESS  OF  CONSTRAINED  LAYER  SYSTEM 

C  ETAS  =  SYSTEM  LOSS  FACTOR 

C  ETA3  =  LOSS  FACTOR  OF  N-TH  CONSTRAINED  LAYER  SYSTEM 

C  ETA3PM  =  LOSS  FACTOR  OF  CONSTRAINING  LAYER  IN  IHE  N-TH  LAYER  (=0 

C  FCP  =  FREQUENCY  OF  THE  DAMPED  PLATE  {HD 

C  IP     =  FREQUENCY  OF  UNDAMPED  TLATE   (HZ) 

C  GC  =  GRAVITATIONAL  CONSTANT 

C  HI  =  THICKNESS  OF  DASE  PLATE 

C  H2  =  THICKNESS  OF  1ST  VISCOELASTIC  LAYER 

C  113  =  THICKNESS  OF  M-TH  CONSTRAINED  LAYER  SYSTEM 

C  IK.  =  TOTAL  DAMPED  PLATE  THICKNESS 

C  II1PRM  =  THICKNESS  OF  BASE  LAYER  IN  N-TH  LAYER 

C  H2PRI1  =  THICKNESS  OF  VEH  III  N-TH  LAYER 

C  H3PRM  =  THICKNESS  OF  CONSTRAINING  LAYER  IN  N-TH  LAYER 

C  KQR  =  HAVE  NUMBER  OF  UNDAMPED  PI  ATE 

C  NU1  =  POISSON'S  RATIO  OF  BASE  PLATE  AND  CONSTRAINING  LAN.ERS 

C  IIU2  =  POISSON'S  RATIO  OF  VISCOELASTIC  MATERIAL 

C  NU2PRM  =  POISSON'S  RATIO  OF  VEM  IN  ll-TH  LAYER 

C  RII01  =  DENSITY  OF  BASE  PLATE  AND  CONSTRAINING  LAYERS 

C  RH02  =  DENSITY  OF  VISCOELASTIC  MATERIAL 

C  RH02PM  =  DENSITY  OF  VEH  IN  N-TH  LAYER 

C  RII03  =  DENSITY  OF  N-IH  LAYER  (  =  RH01) 

C  RH03PM  =  DENSITY  OF  N-TH  LAYER  CONSTRAINING  LAYER 

C  T  =  TEMPERATURE  VARIABLE 

C  UP  =  FREQUENCY  OF  UNDAMPED  PLATE   (RAD/SEC) 

C 

C  C, ETFROL , FROL, FROM, ML , HROM, N , SH , SI , TO  =  COETFICIENIS  FUR  THE 

C  REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS 

C 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

REAL    C,  El  , E3,E3TRM, EHCUBE, ETAS, ET A3, ETA3PM, ET TROL , FCP , FP , FROL 

REAL     FROM,  GC,  HI  ,  H?  ,  II 3  ,  IK. ,  H  1  PRM  ,  II2PRM,  H3PRM,  KQR,  ML  ,  MRUM,  II.NUl 

REAL     NU2,NU2PM,PI  ,  RIIO  1 ,  RII02  ,  RH02PM,  SH  ,  SL  ,  SUD  1  ,  1  ,  T  0  , 1  IP 

REAL  RH03,RH03PM 
C 

DIMENSION  MP( 17 ) 

PI=<i  .OxATAlK  1.0) 
C 

0PEN(NNIT  =  10,FII.  E=  'PLTFRQ',5TATUS='0ID') 

Or  Ell  (UNI  T  =  l  1  ,  r  I  L  E=  '  IIIULYR'  ,  SI  A  I  US  ^  'OLD'  ) 
C 

C  ASSIGN    LAYER     i 'If  I  CKIIF.SSES 
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70 
70 
C 

c 
c 


111  =0.5 

I I  ?  -  U     0  1  5 
II1PRM  =  0  .0937  5 
II2PRM=P  .015 
|I3PRM=0  .  0937  5 

M3  =  II1PRMHI2PRM«M3PRM 

III  =  111  «  H2<  Hi 

ASSIGN    MATERIAL     COIISIANTS 

El =1 . 0E7 
RIIOl  =  0  .  0968 
Hill  =0  .33 
RII02  =  0  .  035 
1(112=0.49 
RMU2PM=0 . 035 
NU2P.I  =  0  .  5 
RIIO3  =  0.  0968 
E3PRII=1 . 0C7 
RII03PM  =  0.  0968 
ETA3f'M=0.0 
ETA3    =0.0 

GC=386  .  0 

DEFIIIE    VISCOELASTIC    COIISTAHTS    FOR    IIOMOGRAN    EQUATIONS 

T  0  =  101.0 

rRon=2 . oei 

MRUM=688  .  91 
tl  =  0.275 
ML  =8.  7 
ETFROLM  .08 
S  L  =  0  .  1  5 
5ll  =  -0.  55 
FROL=5000.0 
C  =  2.5 

ITRITE    PLATE    CHARACTERISTICS 

HRITEC  11  ,700)  '111    =',111 

IIRITEdl  ,700)  '112    =  ',H2 

IIRITEdl, 700)  »  MI  PRIME    =',H1PRH 

IIRI1EC  1  1  ,  700)  'H2PRIME    =',II2PRM 

IIRITEdl, 700)  'H3PRII1E    =',H3PRI1 
IIRI  1E(  1  1  ,  x  J 
IIIUTE(  1  1  ,  x) 

0  IORMATC A12,F6 . 5) 

IIRI  T  EC  11,7  01  )     '  IEMP1  ,  'MODE'  ,  TCP'  ,  'ETAS' 

1  FORMAT* A8, 3X.A5.2A15) 
UR I T  EC  11, X) 

PLA1 E  MODE  LOOP 

DO  100  1=1,17 
READC  10,  x)  IIP <  I  ) 

CALCULATE  MODAL  FREQUENCY  AND  HAVE  NUMBER  0T  UNDAMPED  PIAIE 

SUBl=SQRT((Elx(||<iXX3)xGC)/(12.X(l.-IIUlxx2)xRHOlxH1)) 

KQR  =  IIPd  )/SUDl 
FP  =  IIP(  I)/C2.xri) 

TEMPERATURE  LOOP 

T=30.0 

DO  200  K=l , 15 

COMPUTE    SY51EM    LOSS    FACTOR    AND    STIFFNESS    TOR    N-1H    LAYER 

CUL    RKUCIIIPRM,  M2PRN,  M3PRM,  LI  ,  RHOi  ,  IIU1  ,  IIU2.,i1,  RH02PM,  E3PRI1,  RllUiPfi, 
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CL1  A.',  I'll,  [  ,  [  l\  KQR,  El  AS,  LIICUBC  ,  TCP,  T  0,  FROM,  f  IK  Ml,  II,  ML  ,  El  f  ROL  ,  SL  ,  311, 

c  r  r  o  i  ,  o 

COIIVLRI     RESULTS     FROM    11-111    LAYER    CALCULAIIOII     TO    TOTAL     PLATE 
E1A3  =     0 . 0 

I  i'    =   lire  I )/(  2  .  xpj ) 
f  n  =  n:p 

E3-LIH  UBE-'(  Il3xx3) 

COMPUTE    SYSTEM    LOSS     FACTOR    AMD    TREQUEIICY    TOR    TOTAL     PLATE 

CALL     RKU  (111  ,  112,  H3,  El  ,  RIIOl  ,11111  ,  IIU2,  RII02.  El ,  RII03 ,  ETA3,  T,  PP,  KQR,  ETAS, 
CCHCUDE,  TCP,  T0,FROM,MROM,H,ML  ,  EI  FROL  ,  SL  ,  Sli,  TROL  ,  C  ) 

PRINT  RESULTS 


"RITE(  1  I,  702)  T, I, TCP, ETAS 
7  02    rORMATC  5X,F7.3,2X;I2,3X,2E15.<h) 
C 

C   NEXT  TEMPERATURE 
C 

T  =  T  f  5  .  n 

200       con  r i nut 

c 

C       IICXT     MODE 
C 

HRITEC 1 1, x) 
100  CONTIIIUE 

C 

CLOSEUINI  T  =  l  1  ) 

CLOSE(UII1T  =  10) 

EIID 
C 
C    xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C   XXXXXXXXXXXXXXXXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXXXK*XXXXXXXXXXXXXXXXKXXXX 

C 

SUDROUTKIE    RKU(H1,H2,H3,E1,RH01  , NU1 , NU2, RH02,  E3  ,  RII03,  ETA  3  ,  T  ,  FP  , 
CKQR,ETAS,EHCUBE,  FCP,  TO,  FROM,  MROM,  N,  ML  ,  ETFROL,SL,SM,FROL,C) 
C 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

THIS  SUBROUTINE  CALCULATES  VISCOEIASTIC  PROPERTIES  BASED  ON  THE 
UNIVERSITY  OF  DAYTON  REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS,  AND  THEN 
CALCULATES  PLATE  STIFFNESSES  AND  LOSS  FACTORS  BASED  ON  THE  ROSS- 
KERIIIN-UNGAR  EQUATIONS. 


C 

c 

C 
C 

c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THE  FOLLOWING  ADDITIONAL  VARIABLES  ARE  DEFINED  FOR  USE  IN  THIS 
SUBROUTINE: 

A  =  COEFFICIENT  FOR  NOMOGRAM  EQUATIONS 

ALPHIM, ALPHRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 

ALPHA  IN  THE  RKU  EQUATIONS 
BIM.BRE  =  IMAGINARY  AIID  REAL  COMPONENTS  OF  COEFFICIENT  'B'  IN 

THE  RKU  EQUATIONS 
Ci,D  =  COEFFICIENTS  TOR  RKU  EQUATIONS 
DELIM,DELRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 

DELTA  IN  THE  RKU  EQUA1 IONS 
DEHS  =  COMBINATION  OF  MATERIAL  DENSITIES  USED  TO  COMPUTE  THE 

FREQUEHCY  OF  THE  DAMPED  PLATE  AIID  DAMPING  LAYERS 
EHCUBE  =  EQUIVALENT  STIFFNESS  OF  DAMPED  PLATE  AS  COMPUTED 

USING  RKU  EQUATIONS 
ETA2  =  LOSS  FACTOR  OF  VEM  COMPUTED  IN  NOMOGRAM  EQUATIONS 
ETA210  =  L0G10(ETA2) 

EIAS  =  SYSTEM  LOSS  FACTOR  COMPUTED  BY  RKU  EQUATIONS 
FCP  =  MODAL  FREQUENCY  OF  DAMPED  PLATE  <.\17) 
TP  =  MODAL  FREQUENCY  OF  UNDAMPED  PLAIE  (112) 
TR  =  REDUCED  FREQUENCY  OF  VEM 
FRIO  =  L0G1(1(  FR) 
G2  =  SHEAR  MODULUS  OF  VEM 
1121, H31  =  RKU  EQUATION  COEFFICIENTS 
1110  =  L0G1UCG2)  A3  CUMPUIED  BY  NOMOGRAM  EQUATIONS 


103 


C        l-JCP  =  FREQUENCY  Or  DAMPED  PLAIE  (RAD/SCO 

c 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

REAl     A,  At  rillM,  ALPIIRE,  DIM,BRE,C.C1  ,  D,  DLL  HI,  DEI  RE,  DENS,E1,E2,E3 

REAL  EHCUnE,ETA2,LIA3.F.TA210, CIAS,  E  I  FROL, TCP, FP,rR. FRIO, FRUL,  FROM 

REAL    G.  G2,GC,  Ml  ,112,113,1121  ,M31  ,  KQR,  Ml  0  ,  ML  ,  I1R0M,  II,  NU  1  ,  MU2  ,  P  I  ,  RII01 

REAl     RII02,  RII0  3,5II,SL  ,5Uni  ,  SUBZ,SUD3,  SUD'.,SUD5,  5UC6,  SUD7,  5UD8 

REAL    3UD9,  SUI»10,3Ul)l  1,  T,  1  0,HCP 
C 

ri=<5  .OXATAIK  1.0) 

GC=386 .0 
C 

C   CALCULATE  PROPERTIES  OF  VEI1  FOR  GIVEN  TEMPERATURE  AMD  MODE 
C 
501    rR10  =  LOG10(rP)-(12.x(T-T0))/(525.-»T-T0) 

FR=10.*x(FRl 0) 

A-(f RIO-lOGIOC KROL )  )/C 

SUDl=Cx((SL  +  SM)xA«(3L-5M)K(l.-SQRT(l.-»Axx2)))/2. 

CTA210=LOG10(ETFROL )+SUBl 
C 

C   VISCOELASTIC  LOSS  FACTOR 
C 

ETA2=10.xx(ETA210) 
C 

C   VISCOELASTIC  SHEAR  MODULUS 
C 

Sl)D2  =  2.XLOG10(MROIVHL  ) 

SUB3=1 .+(FROM/FR)XXM 

M10=LOG10(ML)+SUD2/SUD3 

G2=10.xx(M10) 
C 

C   VISCOELASTIC  YOUNG'S  MODULUS 
C 

E2=G2*2.x(l .+NU2) 
C 

C   CALCULATIONS  FOR  DAMPED  PLATE  USING  RKU  EQUATIONS 
C 

G=G2/(E3XH3XH2*KQR) 

H21  =  (Hl-fM2)/2. 

M31=H2+(Hl+H3)/2. 

C1=E1XH1X(1.+G)+GXE3XH3X(1.-ETA2XETA3) 

D  =  GXE1X||]XETA2  +  GXE3XH3X(ETA2  +  ETA3) 

SUD''«  =  GxElxHlxE3xH3x(M3lxx2) 

ALPNRE=SUD^x(Clx(l.-ETA2xETA3)+Dx(ETA2+ETA3)) 

ALPHIM=SUB<ix<ClX(ETA2+ETA3)) 

SUB5  =  ElxmxE2x||2xH31 

BRE=SUB5x(Cl+DxETA2) 

BIM=SUB5X(ClxETAZ-D) 

SUD6=2.XGXE2XH2XE3XH3XH21XM31 

SUD7=l.-2.xETA2xETA3-(ETA2xx2) 

SUna=2.*ETA2+ETA3-(ETA2**2)HETA3 

DELRE=SUB6X(C1XSUB7+DXSUB8) 

DELIM=SUB6X(C1XSUB8-DXSUB7) 
C 

SUB9  =  UZ./(C1XX2  +  DXXZ))X(ALPHRI:-BRE-DELRE) 

EMCUBE=Elx(Mlxx3)«E3x(H3xx3)+SUB9 
C 

C   MODAL  FREQUENCY  OF  DAMPED  PLATE 
C 

DENS=RM0lxHl+RH02xH2+RH03xH3 

SUB10=(EMCUBExGC)/( 1 2 . x ( MU1 xx2 ) xDENS ) 

WCP=KQRXSQRT(SUB10) 

FCP=HCP/(2.XPI ) 

c 

C   COMPARISON  OF  FP  AND  FCP 


C 


IFCABSU  .-FP/FCD  .LE.  0.10)  THEN 

GOIO  500 
ELSE 

FP=FCP 

GUlO  501 
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UNDIF 
COnriME    5Y5TEI1    LOGS    TACIOR 


.00 


5UBll-(12./(ClK*2Ui!<x2))x(ALPHni-DIM-DEHI1) 
ETAS=(l./EHCUnE)*((E3*H3xETA3)»SUbll) 


END 
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APPENDIX  C 

DESIGN  DRAWINGS  FOR  THE  MACHINING 
OF  THE  FLOATING  ELEMENT  AND 
POCKET  PLATE  CONFIGURATIONS 

The  drawings  shown  in  Figures  C.l  and  C.2  were  used  to  machine 
the  pocket  plate  and  floating  element  plate  used  in  the  experiments.  These 
drawings  are  included  to  show  the  relation  of  the  pocket  for  the  ISD— 112  to 
the  cover  plate  and  how  the  viscoelastic  was  protected  from  the  heat  of 
welding. 
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Figure  C.l.   Design  Drawing  of  the  Pocket  Plate  Configuration. 
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Figure  C.2.    Design  Drawing  of  the  Floating  Element  Configuration. 
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APPENDIX  D 

REPRESENTATIVE  MSC/NASTRAN  DATA  DECK  FOR 
THE  DAMPING  CONFIGURATIONS 

This  data  deck  was  used  to  compute  the  modal  frequency  response  of 
the  single  layer  damping  configuration  and  is  a  representative  sample  of 
the  NASTRAN  decks  used  for  the  other  finite  element  models.  The  values 
in  the  damping  table  are  from  a  curve  fit  to  the  modal  loss  factors  estimated 
from  the  modal  strain  energy  method.  Since  the  data  deck  for  the  normal 
mode  and  modal  strain  energy  extraction  is  virtually  identical  to  this  deck, 
the  Case  Control  deck  commands  for  the  normal  mode  analysis  are 
included,  but  are  commented  out. 

The  OUTPUT  request  provides  data  for  an  x-y  plot  of  the  modal 
frequency  response. 

The  units  used  in  this  deck  are  pounds,  inches,  and  seconds. 
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5 


Id  single. ifr 
sol  30 
$ 

5  iittttttitttttttittttttttttittttttttittttiititttttttttitttttittiiitttti 

$      THIS  DECI  IS  TO  COMPUTE  THE  MODAL  FRBQUBNCY  1ESP0ISE  OF  A  15x45 

$  INCH  ALUMINUM  PLATE  WITH  A  CONSTRAINED  VISCOELASTIC  DAMPING  LAKEI. 

$  THE  MODEL  HAS  252  ELEMENTS  HITI  64  ELEMEITS  II  EACH  LAYBfi.  THE  QUAD 

$  ELEMENTS  ABE  OFFSET  FROM  THE  IEI  ELEMEITS  AS  SUGGESTED  BY  THE 

$  LITEBATUIE. 

$ 

$      THE  MATERIAL  PROPERTIES  OF  ISD-112  ARE  FROM  THE  3N  CORPORATION. 

$ 

$      THE  MODAL  LOSS  FACTORS  IN  TBE  DAMPING  TABLE  ARE  FROM  A  CURVB-FIT 

$  TO  TUB  SET  OF  MODAL  LOSS  FACTORS  COMPUTED  FROM  THE  MODAL  STRAIN  ENERGY 

$  METHOD. 

$ 

$      THE  UNITS  USED  IN  THIS  DECK  ABE  POUNDS, INCHES,  AND  SECONDS 

tttttttttttttttttttttittttttttttttttittitttitttttttttttttitttttittttitt 

tlie  3000 

cend 

title  ■-   NODAL  FREQUENCY  RESPONSE  /  200  HZ  /  3M 

let  nod  -  1 

spc  -  1 

dload  =  10 

frequency  -  10 

sdaiplng  -  101 

set  HI  =  95 

svector  -  all 

acceleratlon(plot, phase)  -  HI 

output (iyplot) 

xyprlat  acce  /  95(t3) 

$ 

5  ttiitttttmtn 

$     THE  FOLLOWING  LINES  ARE  THE  CASE  CONTROL  DECI  CARDS  FOR  THE 

$  NORMAL  MODE  EXTRACTION  AND  STRAIN  ENERGY  REQUEST 

$  letbod  -  1 

$  spc  =  1 

$  set  10  -  all 

$  set  11  =  169,tbru,252 

$  ese  *  11 

$ 

$  tttttttitttttit 

BEGIN  BULK 

$  TITLE  =  SINGLE  LAYER  HITH  QUAD  OFFSET 

$  DATA  DECI  PRODUCED  BY  PATNAS  VERSION  2.0:  24-I0V-89  08:31:48 

GRID  1  45.0000  15.0000  0.51500 

GRID  2  45.0000  15.0000  0.50000 

GRID  3  45.0000  12.5000  0.51500 

GRID  4  45.0000  12.5000  0.50000 

GRID  5  45.0000  10.0000  0.51500 

GRID  6  45.0000  10.0000  0.50000 
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GRID 

7 

45.0000  7.50000  0.51500 

GIID 

8 

45.0000  7.50000  0.50000 

GIID 

9 

45.0000  5.00000  0.51500 

GIID 

10 

45.0000  5.00000  0.50000 

G8ID 

11 

45.0000  2.50000  0.51500 

GIID 

12 

45.0000  2.50000  0.50000 

GIID 

13 

45.0000  0.     0.51500 

GIID 

14 

45.0000  0.     0.50000 

GIID 

15 

41.7857  15.0000  0.51500 

GRID 

16 

41.7857  15.0000  0.50000 

GRID 

17 

41.7857  12.5000  0.51500 

GRID 

18 

41.7857  12.5000  0.50000 

GIID 

19 

41.7857  10.0000  0.51500 

GIID 

20 

41.7857  10.0000  0.50000 

GIID 

21 

41.7857  7.50000  0.51500 

GIID 

22 

41.7857  7.50000  0.50000 

GIID 

23 

41.7857  5.00000  0.51500 

GIID 

24 

41.7857  5.00000  0.50000 

GIID 

25 

41.7857  2.50000  0.51500 

GIID 

26 

41.7857  2.50000  0.50000 

GIID 

27 

41.7857  0.     0.51500 

GIID 

28 

41.7857  0.     0.50000 

GIID 

29 

38.5714  15.0000  0.515O0 

GIID 

30 

38.5714  15.0000  0.50000 

GIID 

31 

38.5714  12.5000  0.51500 

GIID 

32 

38.5714  12.5000  0.50000 

GIID 

33 

38.5714  10.0000  0.51500 

GIID 

34 

38.5714  10.0000  0.50000 

GIID 

35 

38.5714  7.50000  0.51500 

GIID 

36 

38.5714  7.50000  0.50000 

GIID 

37 

38.5714  5.00000  0.51500 

GIID 

38 

38.5714  5.00000  0.50000 

GIID 

39 

38.5714  2.50000  0.51500 

GIID 

40 

38.5714  2.50000  0.50000 

GIID 

41 

38.5714  0.     0.51500 

GIID 

42 

38.5714  0.     0.50000 

GIID 

43 

35.3571  15.0000  0.51500 

GIID 

44 

35.3571  15.0000  0.50000 

GIID 

45 

35.3571  12.5000  0.51500 

GIID 

46 

35.3571  12.5000  0.50000 

GIID 

47 

35.3571  10.0000  0.51500 

GIID 

46 

35.3571  10.0000  0.50000 

GIID 

49 

35.3571  7.50000  0.51500 

GIID 

50 

35.3571  7.50000  0.50000 

GIID 

51 

35.3571  5.00000  0.51500 

GIID 

52 

35.3571  5.00000  0.50000 

GIID 

53 

35.3571  2.50000  0.51500 

GIID 

54 

35.3571  2.50000  0.50000 

GIID 

55 

35.3571  0.     0.51500 

GIID 

56 

35.3571  0.     0.50000 

GIID 

57 

32.1429  15.0000  0.51500 
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GRID 

58 

GIID 

59 

GIID 

60 

GIID 

61 

GRID 

62 

GIID 

63 

GRID 

64 

GRID 

65 

GRID 

66 

GRID 

67 

GRID 

68 

GRID 

69 

GRID 

70 

GRID 

71 

GRID 

72 

GRID 

73 

GRID 

74 

GRID 

75 

GRID 

76 

GRID 

77 

GRID 

78 

GRID 

79 

GRID 

80 

GRID 

81 

GRID 

82 

GRID 

83 

GRID 

84 

GRID 

85 

GRID 

86 

GRID 

87 

GRID 

88 

GRID 

89 

GRID 

90 

GRID 

91 

GRID 

92 

GRID 

93 

GRID 

94 

GRID 

95 

GRID 

96 

GRID 

97 

GRID 

98 

GRID 

99 

GRID 

100 

GRID 

101 

GRID 

102 

GRID 

103 

GRID 

104 

GRID 

105 

GRID 

106 

GRID 

107 

GRID 

108 

32.1429  15.0000  0.50000 
32.1429  12.5000  0.51500 
32.1429  12.5000  0.50000 
32.1429  10.0000  0.51500 
32.1429  10.0000  0.50000 
32.1429  7.50000  0.51500 
32.1429  7.50000  0.50000 
32.1429  5.00000  0.51500 
32.1429  5.00000  0.50000 
32.1429  2.50000  0.51500 
32.1429  2.50000  0.50000 
32.1429  0.  0.51500 
32.1429  0.  0.50000 
28.9286  15.0000  0.51500 
28.9286  15.0000  0.50000 
26.9286  12.5000  0.51500 
28.9286  12.5000  0.50000 
28.9286  10.0000  0.51500 

28.9286  10.0000  0.50000 
28.9286  7.50000  0.51500 
28.9286  7.50000  0.50000 
28.9286  5.00000  0.51500 
28.9286  5.00000  0.50000 
28.9286  2.50000  0.51500 
28.9286  2.50000  0.50000 
28.9286  0.  0.51500 
28.9286  0.  0.50000 
25.7143  15.0000  0.51500 
25.7143  15.0000  0.50000 
25.7143  12.5000  0.51500 
25.7143  12.5000  0.50000 
25.7143  10.0000  0.51500 
25.7143  10.0000  0.50000 
25.7143  7.50000  0.51500 
25.7143  7.50000  0.50000 
25.7143  5.00000  0.51500 
25.7143  5.00000  0.50000 
25.7143  2.50000  0.51500 
25.7143  2.50000  0.50000 
25.7143  0.  0.51500 
25.7143  0.  0.50000 
22.5000  15.0000  0.51500 
22.5000  15.0000  0.50000 
22.5000  12.5000  0.51500 
22.5000  12.5000  0.50000 
22.5000  10.0000  0.51500 
22.5000  10.0000  0.50000 
22.5000  7.50000  0.51500 
22.5000  7.50000  0.50000 
22.5000  5.00000  0.51500 
22.5000  5.00000  0.50000 


112 


G8ID 

109 

G8ID 

110 

Glib 

111 

GBID 

112 

GklD 

113 

CUD 

114 

GRID 

115 

GIID 

116 

GBID 

117 

GIID 

118 

GBID 

119 

GIID 

120 

GIID 

121 

GIID 

122 

GIID 

123 

GIID 

124 

GIID 

125 

GIID 

126 

GIID 

127 

GIID 

128 

GIID 

129 

GIID 

130 

GBID 

131 

GIID 

132 

GIID 

133 

GBID 

134 

GIID 

135 

GIID 

136 

GIID 

137 

GIID 

138 

GIID 

139 

GIID 

140 

GIID 

141 

GIID 

142 

GIID 

143 

GBID 

144 

GBID 

145 

GIID 

146 

GIID 

147 

GIID 

148 

GIID 

149 

GIID 

150 

GIID 

151 

GIID 

152 

GIID 

153 

GIID 

154 

GIID 

155 

GIID 

156 

GIID 

157 

GIID 

156 

GIID 

159 

22.5000  2.50000  0.51500 
22.5000  2.50000  0.50000 
22.5000  0.  0.51500 
22.5000  0.  0.50000 
19.2857  15.0000  0.51500 
19.2857  15.0000  0.50000 
19.2857  12.5000  0.51500 
19.2857  12.5000  0.50000 
19.2857  10.0000  0.51500 
19.2857  10.0000  0.50000 
19.2857  7.50000  0.51500 
19.2857  7.50000  0.50000 
19.2857  5.00000  0.51500 
19.2857  5.00000  0.50000 
19.2857  2.50000  0.51500 
19.2857  2.50000  0.50000 
19.2857  0.  0.51500 
19.2857  0.  0.50000 
16.0714  15.0000  0.51500 
16.0714  15.0000  0.50000 
16.0714  12.5000  0.51500 
16.0714  12.5000  0.50000 
16.0714  10.0000  0.51500 
16.0714  10.0000  0.50000 
16.0714  7.50000  0.51500 
16.0714  7.50000  0.50000 
16.0714  5.00000  0.51500 
16.0714  5.00000  0.50000 
16.0714  2.50000  0.51500 
16.0714  2.50000  0.50000 
16.0714  0.  0.51500 
16.0714  0.  0.50000 
12.8571  15.0000  0.51500 
12.8571  15.0000  0.50000 
12.8571  12.5000  0.51500 
12.8571  12.5000  0.50000 
12.8571  10.0000  0.51500 
12.8571  10.0000  0.50000 
12.8571  7.50000  0.51500 
12.8571  7.50000  0.50000 
12.8571  5.00000  0.51500 
12.8571  5.00000  0.50000 
12.8571  2.50000  0.51500 
12.8571  2.50000  0.50000 
12.8571  0.  0.51500 
12.8571  0.  0.50000 
9.64286  15.0000  0.51500 
9.64286  15.0000  0.50000 
6.42857  15.0000  0.51500 
6.42657  15.0000  0.50000 
3.21429  15.0000  0.51500 
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GIID  160  3.21429  16.0000  0.50000 

GIID  161  0.     15.0000  0.51500 

GfilD  162  0.     15.0000  0.50000 

GIID  163  9.64286  12.5000  0.51500 

GEID  164  9.64286  12.5000  0.50000 

GIID  165  9.64286  10.0000  0.51500 

GIID  166  9.64286  10.0000  0.50000 

GIID  167  9.64286  7.50000  0.51500 

GIID  168  9.64286  7.50000  0.50000 

GIID  169  9.64286  5.00000  0.51500 

GIID  170  9.64286  5.00000  0.50000 

GIID  171  9.64286  2.50000  0.51500 

GIID  172  9.64286  2.50000  0.50000 

GIID  173  9.64286  0.     0.51500 

GIID  174  9.64286  0.     0.50000 

GIID  175  6.42857  12.5000  0.51500 

GIID  176  6.42857  12.5000  0.50000 

GIID  177  3.21429  12.5000  0.51500 

GIID  178  3.21429  12.5000  0.50000 

GIID  179  0.     12.5000  0.51500 

GIID  180  0.     12.5000  0.50000 

GIID  161  6.42857  10.0000  0.51500 

GIID  182  6.42857  10.0000  0.50000 

GIID  183  6.42857  7.50000  0.51500 

GIID  184  6.42857  7.50000  0.50000 

GIID  185  6.42857  5.00000  0.51500 

GIID  186  6.42857  2.50000  0.51500 

GIID  187  6.42857  5.00000  0.50000 

GIID  188  6.42857  2.50000  0.50000 

GIID  189  6.42857  0.     0.51500 

GIID  190  6.42857  0.     0.50000 

GIID  191  3.21429  10.0000  0.51500 

GIID  192  3.21429  10.0000  0.50000 

GIID  193  0.     10.0000  0.51500 

GIID  194  0.     10.0000  0.50000 

GIID  195  3.21429  0.     0.51500 

GIID  196  3.21429  0.     0.50000 

GIID  197  0.     0.     0.51500 

GEID  198  0.     0.     0.50000 

GIID  199  3.21429  7.50000  0.51500 

GIID  200  3.21429  5.00000  0.51500 

GIID  201  3.21429  2.50000  0.51500 

GIID  202  3.21429  7.50000  0.50000 

GIID  203  3.21429  5.00000  0.50000 

GIID  204  3.21429  2.50000  0.50000 

GIID  205  0.     7.50000  0.51500 

GIID  206  0.     5.00000  0.51500 

GIID  207  0.     2.50000  0.51500 

GIID  208  0.     7.50000  0.50000 

GIID  209  0.    2.50000  0.50000 

GIID  210  0.     5.00000  0.50000 
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